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and  at  depth  is  important  because  near  surface  cur¬ 
rents  are  “contaminated”  by  local  surface  winds 
which  make  a  determination  of  the  mean  flow  dif¬ 
ficult. 

The  dual-frequency  microwave  current  meter 
has  been  successfully  tested  from  towers  in  the  Atlan¬ 
tic  Ocean  and  in  the  North  Sea.  System  tests  using  a 
mobile  platform  (a  research  vessel  the  AGOR-16 
USNS  Hayes)  are  planned  for  this  July  in  both  the 
Gulf  Stream  and  in  the  shoals  near  Nantucket  Island. 
■ 

(Dale  L.  Schuler,  NRL) 

Wiiiter*Season  Thermal  Stress  to  High- 
Latitude  Reef  Systems 

The  effects  of  cold-air  outbreaks  over  high- 
latitude  reef  systems  are  being  investigated  by 
scientists  at  the  Coastal  Studies  Institute  of  Louisiana 
State  University.  Nan  Walker,  funded  by  a  NASA 
Graduate  Student  Research  Grant,  and  Drs.  H.  H, 
Roberts,  L.  J.  Rouse,  Jr.,  and  O.K.  Huh,  funded  by 
ONR,  have  found  that  these  outbreaks  over  the  reefs 
and  associated  shallow  environments  can  create 
thermal  stress  conditions  severe  enough  to  kill  reef¬ 
building  corals,  thereby  affecting  long-term  reef 
development.  Over  the  last  five  years  digital  thermal 
infrared  data  from  high-resolution  radiometers  of 
the  NOAA-5  and  NOAA-6  environmental  satellites 
have  provided  a  means  for  determining  regional  sea- 
surface  temperature  variations  during  particularly 
cold  outbreaks  of  polar  air.  Synoptic  satellite  data 
were  used  in  conjunction  with  a  computerized  heat 
flux  model  to  investigate  shallow-water  chilling 
processes.  The  model  incorporates  hourly 
meteorological  parameters  and  the  bulk  aerodynamic 
equations  to  hindcast  changes  in  water  temperature 
throughout  critical  cold-air  outbreaks.  These  results, 
coupled  with  in  situ  temperatures  and  biological 
observations  (from  Florida  Bay  and  the  northern 
Bahama  Banks),  have  provided  a  data  base  from 
which  a  comprehensive  appraisal  of  thermal  stress 
events  on  reef  communities  has  been  formulated. 

During  January  1977  and  1981  consecutive  cold 
fronts  crossed  south  Florida  and  the  nt  them 
Bahamas,  causing  abnormal  chilling.  Field  and 
laboratory  experiments  identify  16°C  as  a  thermal 
stress  threshold  for  most  reef-building  corals,  below 
which  they  have  difficulty  maintaining  life¬ 
supporting  processes.  Water  temperatures  below 
16'C  were  experienced  in  Florida  Bay  for  8  con¬ 
secutive  days  during  the  1977  event  and  for  10  days  in 


1981.  Minimum  Florida  Bay  temperatures  obtained 
from  corrected  satellite  data  and  confirmed  from 
“surface  truth”  measurements  were  12.9°C  and 
9.0®C,  respectively.  The  figure  above  illustrates  a 
profile  of  radiation  temperatures  extracted  from 
digital  satellite  data  along  the  transect  A-A’. 

Massive  fish  kills  and  numerous  coral  mor¬ 
talities  were  observed  in  the  study  area  during  the 
1977  and  1981  cold-air  outbreaks.  Results  show  that 
cold-water  stress  conditions  can  exist  over  vast 
shallow-water  areas  with  residence  times  of  several 
days.  Offbank  movement  of  chilled  waters,  driven  by 
northerly  winds,  density  gradients,  and  tidal  pum¬ 
ping,  can  limit  coral  growth  and  control  the  spatial 
distribution  of  reef  development  in  high  latitudes.  ■ 


The  COj  laser  is  focused  on  the  water  surface 
causing  an  explosion  called  “thermal  blow-off.” 
Preliminary  st\>dies  have  indicated  that  the  explosion 
creates  an  acoustic  field  in  the  water  having  a  pulse 
width  of  20-30  ;«ecs.  and  a  frequency  spectrum 
peaking  between  25-30  kHs.  On  page  S  is  an  example 
of  a  typical  data  return  from  one  pulse  of  the  radar. 
The  air  blast  (surface  explosion  noise)  and  the 
bottom  bounce  signals  are  quite  evident.  The  depth  is 
easily  calculated  because  the  difference  between  the 
speed  of  sound  in  water  and  air  is  quite  exaggerated. 
The  measurement  was  taken  while  at  sea  and  un¬ 
derway.  ■ 


(Paul  Twitchell,  ONR  Boston) 
Con  tin  ued  from  page  5 
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(James  S.  Bailey,  ONR) 
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Remote  Sensing  Research 

by  Paul  F.  Twltchall 
Offlc*  of  Naval  ftotoarch,  Bottoo 


The  U.S.  Navy  has  a  long  history  in  remote  sens¬ 
ing  research  and  its  application  to  communica¬ 
tion,  navigation,  surveillance,  weapons  detection  and 
delivery  systems.  The  Geography  Division  of  the  Of¬ 
fice  of  Naval  Research  (ONR)  initiated  in  the  1950’s 
a  research  program  focusing  on  remotely  sensing  the 
environment.  ONR  was  a  pioneer  in  this  field  and 
continues  today  as  an  important  force  in  many  new 
areas  of  remote  sensing  research.  In  1962  ONR  spon¬ 
sored  the  first  International  Symposium  on  the 
Remote  Sensing  of  the  Environment.  In  recent  years 
these  meetings  have  grown  into  major  conferences 
each  resulting  in  proceedings  of  several  volumes.  In 
1973  ONR’s  Remote  Sensing  Program  was  reviewed 
in  Naval  Research  Reviews'.  In  this  issue  are  four 
papers  focusing  on  a  few  facets  of  remote  sensing  of 
the  ocean. 

For  more  than  forty  years,  Navy  scientists  have 
been  integrating  advances  in  microwave  technology 
into  the  fleet  including  shipboard  and  aircraft  fire 
control  radars,  weapons  guidance,  and  ocean  surveil¬ 
lance  from  space.  In  the  past  twenty  years  the  use  of 
radio  waves  to  measure  ocean  surface  dynamics  has 
contributed  greatly  to  the  understanding  of  wave 
dynamics.  Navy  scientists,  particularly  the  late  John 
W.  Wright  of  the  Naval  Research  Laboratory,  have 
played  a  significant  rok  in  studying  phenomena  at 
the  air-sea  interface  relating  to  short-wave  modula¬ 
tion  by  long  waves,  the  spectral  development  of 
ocean  surface  waves,  and  wave  growth  predictions. 
Many  of  the  basic  scientific  questions,  however,  have 
not  been  answered;  for  example,  phenomena  such  as 
nonlinear  wave-wave  energy  transfer  and  strong  non¬ 
linear  wave-current  interaction  require  improved 
theoretical  understanding  to  fully  explain  radar 
amplitude  modulation  and  to  extract  wave  spectra  in¬ 
formation  from  radar  returns.  Professor  Erik  Mollo- 
Christensen’s  article  in  this  issue  addresses  some  of 
these  unanswered  questions. 

Microwave  remote  sensing  of  the  ocean  uses 
radiometric  and  altimetry  devices.  For  about  a 
decade  scanning  microwave  radiometers  have  flown 
on  U.S.  satellites.  The  Navy  has  and  continues  to 


develop  techniques  for  using  the  "brightness  temper¬ 
ature”  as  measured  by  the  radiometers.  One  practical 
application  has  been  the  identification  and  classifica¬ 
tion  of  sea  ice.  Among  the  information  gained  from 
radar  altimeters  on  satellites  are  the  locations  of  large 
ocean  eddies  and  the  slope  of  the  ocean  surface 
across  a  major  current  system.  The  location  and 
speed  of  currents  can  be  estimated  from  space  using 
ocean  slope  measurements. 

The  infrared  sensors  that  are  common  today  on 
aircraft  and  satellites  are  traceable  to  basic  research 
funded  in  part  by  the  physics  and  chemistry  branches 
of  the  Office  of  Naval  Research  in  the  1950’s.  Hard¬ 
ware  developments  and  their  applications  were  ac¬ 
celerated  through  the  Infrared  Information  Sym¬ 
posiums  organized  and  run  by  ONR.  Special  credit 
should  be  given  to  two  late  ONR  scientists,  T.  B. 
Dowd  of  the  Boston  Branch  Office  and  L.  B.  White 
of  the  Chicago  Branch  Office  for  directing  these  sym¬ 
posia  which  provided  a  forum  for  the  development  of 
modern  electrooptical  weapons  systems. 

H.  L.  Clark  and  other  scientists  at  the  Naval  Re¬ 
search  Laboratory  developed  an  early  (1942)  infrared 
system  for  remotely  sensing  the  sea  surface  from 
lighter-than-air  shipsb  By  the  1960’s  airborne  in¬ 
frared  systems  were  in  common  use.  Two  ONR 
funded  investigators,  P.  M.  Saunders  at  Woods  Hole 
Oceanographic  Institution  and  E.  D.  McAllister  at 
Scripps  Institute  of  Oceanography,  used  airborne  in¬ 
frared  systems  in  their  basic  investigations  of  oceanic 
processes.  Contemporary  ONR-funded  efforts  that 
employ  infrared  sensors  are  reported  in  this  issue  by 
Otis  B.  Brown  and  Robert  Evans  of  the  University  of 
Miami  and  Oscar  Huh  and  Lawrence  Rouse  of  Loui¬ 
siana  State  University. 

By  using  environmental  satellites  in  the  early 
1960’s,  the  Navy  pioneered  in  research  which  ob¬ 
tained  on  a  synoptic  (l(X)0’s  km)  and  meso  (lOO’s  km) 
scale  environmental  characteristics  of  the  ocean  and 
atmosphere.  Over  the  past  two  decades  advances  in 
quantifying  satellite  data  have  resulted  in  cost  saving 
for  data  acquisition,  and  a  better  understanding  of 
the  physical  processes  of  the  ocean.  The  Office  of 
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Naval  Research  has  sponsored  many  scientific  meet¬ 
ings  and  workshops  to  advance  the  application  of 
remote  sensing  tools  for  the  study  of  oceanic  proc¬ 
esses.  At  the  meeting  on  Physical  Oceanography  and 
Satellites  held  at  Wooclshole  in  June  1978,  Professor 
Verner  Suomi  concluded  with  the  following 
statement’: 

Remote  sensing  and  data  collection  from 
satellite  platforms  show  considerable  promise 
for  improving  our  understanding  of  a  wide 
range  of  scales  of  ocean  flows  and  their  in¬ 
teraction  with  the  atmosphere  and  underlying 
surface.  The  data  volume  from  these  new 
tools  will  be  immense  compared  to  traditional 
methods  of  gathering  oceanographic  data.  If 
the  full  benefits  of  these  new  tools  and  data 
sets  are  to  be  realized,  new  methods  are  re¬ 
quired  to  convert  these  composite  data  sets 
into  information.  The  problem  is  similar  to 
that  experienced  by  the  meteorological 
satellite  operations.  We  therefore  recommend 
that  rapid  analysis  tools  using  mintcomputers 
and  interactive  video  displays  be  made  avail¬ 
able  to  oceanographic  researchers:  further 
that  a  modest  number  of  these  tools  be  made 
available  to  individual  researchers  or  to  small 
groups  rather  than  through  a  simple  or  small 
number  of  large  institutional  data  processing 
facilities.  Our  present  ability  to  fully  utilize 
these  vast  data  sets  is  limited.  Nevertheless,  we 
recommend  that  as  many,  and  as  complete, 
data  sets  as  possible  be  archived,  taking  into 
account  reasonable  cost  but  at  the  same  time 
being  mindful  of  the  future  potential  value  of 
the  data  for  studies  of  the  time-history  of  large 
scale  flows.  We  stress  that  any  archiving 
scheme  will  ha  ..lue  only  in  direct  propor¬ 
tion  to  provisions  for  economical  and  ready 
access  to  the  data  of  the  research  and  opera¬ 
tional  community. 

The  Navy  remote  sensing  program  is  addressing 
the  needs  suggested  by  Professor  Suomi  and  Navy 
environmental  support  requirements  outlined  in  the 
following  table. 

As  indicated  by  Table  1  the  Navy  environmental 
support  requirements  span  a  broad  range  of  time-and 
space-scales  demanded  by  each  of  the  several  Navy 
mission  areas.  Applications  of  satellite  oceano¬ 
graphic  data  include  support  of  fleet  mission  oper¬ 
ations  and  support  of  various  research  and  develop¬ 
ment  programs. 

There  exists  a  Navy  Operational  Requirement 


Table  I 

Navy  Environmental  Support  Requiremenls 
(ADAPTED  FROM  Noble  and  Fell) 


MISSION  SUPPORT: 

STRIKE,  ANTISUBMARINE, 
ANTISURFACE,  AMPHIBIOUS,  and 
MINE  WARFARE 

TIME-SCALES 

HOURS  (AIR  OPERATIONS) 

DAYS  (TYPHOON  FORECASTING, 

ICE  OPERATIONS) 

WEEKS  (MISSION  PLANNING,  SHIP 

ROUTING) 

SPACE-SCALES 

10m  (COASTAL  PROCESSES) 

1-50km  (OCEANIC  FRONTS,  EDDIES, 

CURRENTS) 

50-1000km  (SEVERE  STORMS,  GLOBAL 

“WEATHER”) 

APPLICATIONS 

TACTICAL  AND  STRATEGIC  MISSION 
OPERATIONS 

RESEARCH  AND  DEVELOPMENT 
PROGRAMS 


for  Satellite  Measurement  of  Oceanographic  Param¬ 
eters  (SMOP  0R-W0527-0S).  The  specifications  for 
horizontal  resolution,  precision,  and  accuracy  of 
measurement  of  identified  parameters  are  shown  in 
Table  2.  The  dual  figures  in  the  statement  of  reso¬ 
lution,  precision,  and  accuracy  represent  estimates  of 
minimal  acceptable  thresholds  for  operational  utility 
when  the  document  was  prepared*. 

To  achieve  the  operational  requirements  out¬ 
lined  above  there  is  a  need  for  basic  and  applied 
research  on  sensors,  image  enhancement  and  display, 
data  interpretation,  integration  of  satellite  data  into 
numerical  models  and  other  topics.  The  paper  in  this 
issue  by  Thomas  Brubecker  of  Colorado  State  Uni¬ 
versity  is  an  example  of  the  efforts  to  translate  in 
near  real  time  the  immense  volumes  of  data  stream¬ 
ing  down  from  satellites  into  a  useful  presentation 
for  an  operational  person  such  as  an  aircraft  pilot.  In 
a  later  issue  of  Naval  Research  Reviews  there  will  be 
additional  articles  on  improving  satellite  data  trans¬ 
mission  and  processing. 

Another  area  of  research  is  the  application  of 
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Table  2 

Operalionul  Requiremenls  Sulellile  Measurements  of  Oceanographic  Parameters 
ISMOP  OP-  W0527-0Si 


PARAMETER 

HORIZONTAL 

RESOLUTION* 

PRECISION 

ACCURACY* 

SURFACE  WIND 

SPEED 

DIRECTION 

10/ 125  km 
10/ /50  km 

5%l  120% 
±5°/  /lO* 

2m/sec/  /4m/sec 
±10°/  /22.5° 

SEA  SURFACE 

TEMPERATURE 

10/ 125  km 

0.25/  /0.8°C 

0.5/  /1.0°C 

SEA  VERTICAL  TEMP. 
STRUCTURE** 

10/ /1 00  km 

0.25/  /1.0°C 

1.0/  /2.0°C 

WAVES 

SIGNIFICANT  HEIGHT 
AMPLITUDE  COMPONENTS 
WAVELENGTH  COMPONENTS 
DIRECTION 

10/ 125  km 

1/  /25  km 

1/  /25  km 

1/  /50  km 

0.3m/  /10% 
0.3m/  /0.7m 
5%/  /1 5% 

10°/  /30° 

0.3  ml  n0% 

0.3m/  /0.7m 

5%/ 135% 

10/  /45° 

ICE 

COVER 

THICKNESS 

AGE 

0.5/ 125  km 

21 150  km 
10/  /50  km 

10%/  /25% 
0.25/  /2m 

6/  /1 2  mos. 

12%/ 130% 

0.51 12m 

61 132  mos. 

ICEBERGS 

0.015/  /0.1  km 

N/A 

POSITION  0.5/ 12  km 

•  OPERATIONAL  TARGET  / 1  MINIMUM  ACCEPTABLE  THRESHOLD 
•*  DEPTH  300/ /50m;  VERTICAL  RESOLUTION  ±2/ /10m 

satellite  imagery  of  ocean  surface  characteristics  and 
atmospheric  profiling  into  numerical  models  of 
oceanic  and  marine  boundary  layer  dynamics.  At  the 
present  time  there  is  a  need  to  find  resources  to 
continue  and  expand  efforts  to  optimize  the  ingestion 
of  a  sateiKte  data  into  numerical  models  used  to 
describe  and  predict  operationally  important  envi¬ 
ronmental  information. 

The  remote  sensing  articles  in  this  issue  represent 
only  a  small  fraction  of  the  ongoing  remote  sensing 
programs  of  the  Navy.  In-house  remote  sensing  ac¬ 
tivities  of  the  Navy  include  research  at  Naval  Envi¬ 
ronmental  Prediction  Research  Facility  (NEPRF) 
where  significant  contributions  have  been  made  on 
tactical  applications  of  satellite  data  and  the  develop¬ 
ment  of  a  Satellite  Processing  and  Display  System 
(SPADS).  At  the  Naval  Research  Laboratory 
pioneering  work  on  sensor  technology  is  interna¬ 
tionally  recognized.  The  Naval  Surface  Weapon’s 


Center,  Dahlgren  (NSWC/DL)  and  the  Navy  Ocean¬ 
ographic  Office  (NAVOCEANO)  are  pursuing 
research  using  satellite  data.  At  the  Naval  Ocean 
Research  and  Development  Activity  (NORDA)  scien¬ 
tists  are  developing  techniques  for  effective  opera¬ 
tional  use  of  satellite  imagery  and  investigating 
methods  for  applying  satellite  data  to  numerical 
models.  Other  in-house  efforts  related  to  systems 
developments  include  the  multispectral  satellite  data 
analysis  at  Naval  Ocean  Systems  Center  (NOSC)  and 
the  application  of  satellite  sea  surface  temperature 
data  at  Naval  Underwater  Systems  Center  (NUSC). 
Scientists  of  the  West  Coast  ONR  Office  have  a  long 
tradition  in  fostering  new  remote  sensing  cooperative 
efforts  with  other  agencies  and  various  elements  of 
the  operational  Navy. 

Navy  scientists  and  their  university  colleagues 
funded  by  the  Office  of  Naval  Research  and  Naval 
Air  Systems  Command  (Air  370)  Atmospheric/ 
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Oceanographic  Program  communicate  with  each 
other  through  professional  societies.  On  occasion 
workshops  or  meetings  are  organized  by  ONR  to  en¬ 
courage  advances  in  areas  of  research  of  interest  to 
the  Navy.  The  results  of  ONR-sponsored  activities 
are  well  promulgated  through  professional  society 
publications  or  special  volumes:  recent  examples  in¬ 
clude  a  workshop  report  in  the  Bulletin  of  the  Ameri¬ 
can  Meteorological  Society’  and  a  book*.  Navy  plan¬ 
ners  and  managers  draw  upon  such  resources 
through  groups  such  as  the  Navy  Environmental 
Remote  Sensing  Coordinating  Advisory  Committee 
(NERSCAC)  to  introduce  research  products.  For  ex¬ 
ample,  a  specific  research  finding  could  be  integrated 
into  the  Tactical  Environmental  Support  System 
(TESS)  development.  The  Navy  operational  com¬ 
munity  is  presently  formulating  plans  for  effective 
and  affordable  impiementation  of  satellite  products 
into  global  and  regional  environmental  prediction 
systems. 

The  author  wishes  to  acknowledge  contributions 
from  colleagues  within  the  Office  of  Naval  Research, 
the  Naval  Research  Laboratory,  and  from  other 
organizations  mentioned  throughout  the  text.  ■ 
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New  Laser  Technique  for  Inshore  Bathj^metri^ 


Rapid  surveillance  of  inshore  water  for  bath¬ 
ymetry  and  underwater  obstacles  is  important 
for  a  variety  of  applications.  Extensive  research 
has  been  performed  in  this  area  during  the  past 
several  years  using  a  pulsed  blue-green  laser  to 
map  shallow  water  depth  from  aircraft. 
Limitations  to  this  system  such  as  turbidity  have 
necessitated  research  on  a  more  versatile  or  less 
restricting  system. 

A  new  laser  technique  is  being  studied  for  in¬ 
shore  bathymetry  and  object  detection.  The 
system  is  a  hybrid  system  consisting  of  an  infrared 
laser  (CO  ’  )  transmitter  and  a  highly  sensitive 
microphone  receiver  that  is  located  in  the  air. 

Continued  on  page  6! 
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ACRONYMS 

These  Acronyms  Appear  Frequently  in  This  Issue 

AVHRR 

Advanced  Very  High  Resolution  Radiometer 

CNOC 

Commander  Naval  Ocaanography  Command 

CZCS 

Coastal  Zone  Color  Scanner 

CMSP 

Defense  Meteorological  Satellite  Program 

ERTS 

Earth  Resources  Technology  Satellite 

FOV 

Field  of  View 

GOSSTCOMP 

Global  Operational  Sea  Surface  Temperature  Computation 

HIRS 

HIgh-Resolutlon  Infrared  Sounder 

IFOV 

Instantaneous  Field  of  View 

INDEX 

Indian  Ocean  Experiment 

IR 

Infrared 

IRIS 

Infrared  Interferometer  Spectrometer  (also  Infrared  Information 

LANOSAT 

Symposium) 

Land  Satellite  (Earth  Resources  Technology  Satellite) 

MclDAS 

Man-computer  Interactive  Data  Access  System 

NASA 

National  Aeronautics  and  Space  Administration 

NAVAIR 

Naval  Air  Systems 

NAVOCEANO 

Navy  Oceanographic  Office 

NEPRF 

Naval  Environmental  Prediction  Research  Facility 

NERSCAC 

Navy  Environmental  Remote  Sensing  Coordinating  Advisory 

NESS 

Committee 

National  Earth  Satellite  Service 

NIMBUS 

Polarorbiting  Meteorological  Satellite 

NOAA 

National  Oceanic  and  Atmospheric  Administration 

NOAA-5,  NOAA-6,  NOAA-7,  «tc. 

Operational  Polarorbiting  Meteoroiogical  Satellites 

NOMSS 

Naval  Oceanographic  and  Meteorological  Support  System 

NOROA 

Naval  Oceanographic  Research  and  Development  Activity 

NOSC 

Naval  Ocean  Systems  Center 

NRL 

Naval  Research  Laboratory 

NSWC 

Naval  Surface  Weapons  Center 

NVSC 

Naval  Underwater  Systems  Center 

RSMAS 

Rosentlel  School  of  Marine  and  Atmospheric  Science 

SAR 

Synthetic  Aperature  Radar 

SASS 

SEASAT-A  Scatterometer  System 

SCAMS 

Scanning  Microwave  Radiometer 

SEASAT 

Sea  Satellite  (a  satellite  designed  to  measure  oceanic 

SEBEX 

characteristics,  it  operated  for  a  short  time  in  1978) 

Surface  Expression  of  Bathemetry  Experiment 

SLAR 

Side  Looking  Airborne  Radar 

SMMR 

Scanning  Multifrequency  Microwave  Radiometer 

SST 

Sea  Surface  Temperature 

TESS 

Tactical  Environmental  Support  System 

THIR 

Temperature,  Humidity,  Infrared  Radiometer 

VAS 

Vertical  Atmospheric  Sounder 

VTPR 

Vertical  Temperature  Profile  Radiometer 

Visible  and  Infrared 
Satellite  Remote  Sensing; 
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greater  importance  to  the  understanding  of  oceanic- 
variability  than  the  quasi-quantitative  values  of  the 
data. 

Satellite-based  remote  sensing  of  the  environ¬ 
ment  first  reached  public  notice  with  display  of  color 
earth  pictures  taken  from  the  Mercury  manned 
spacecraft.  In  fact  scientific  sounding  rockets  and 
military  passive  sensors  were  routinely  mapping  the 
oceans  before  this  time'  ■’ '.  Multi-spectral  color  pic¬ 
tures  captured  public  and  scientific  fancy.  However, 
it  was  to  be  some  15  years  later  with  the  launch  of  the 
NIMBUS-7  Coastal  Zone  Color  Scanner  before  these 
promises  could  start  to  be  realized  over  oceanic 
areas;  the  Earth  Resources  Technology  Satellite,  later 
called  LANDSAT,  provided  quasi-operational  land 
remote  sensing  in  the  early  1 970s. 

Application  of  infrared  remote  sensing  to  the 
study  of  meteorological  variability  was  quickly  ap¬ 
preciated  by  analysts  and  an  operational  military 
program  was  initiated  in  1963.  This  program  was 
managed  by  the  U.S.  Air  Force  and  concentrated  on 
determination  of  world-wide  cloudiness,  cloud-top 
temperatures  and  tactical  applications  for  aircraft 
logistical  use.  It  was  not  focussed  on  the  study  of 
oceanic  surface  thermal  variability.  The  past  decade 
has  seen  active  Navy  participation  in  these  efforts, 
now  called  Defense  Meteorological  Satellite  Pro¬ 
gram,  with  a  gradual  orientation  of  the  system 
toward  improving  ocean  thermal-imaging  capabil¬ 
ities. 

In  this  report  we  will  discuss  the  visible/infrared 
sensing  problem,  some  applications  and  conclusions. 
The  sensing  problem  logically  splits  into  two  sec¬ 
tions:  technology  to  acquire  the  data,  ie.,  sensors, 
and  algorithms  or  processing  paths  for  the  data  once 
acquired.  We  will  not  discuss  the  sensor  problem  ex¬ 
cept  as  needed  for  algorithm  development. 

Sensor  Processing 

Sensor  data  processing  can  be  divided  into  the 
following  4  categories; 

1)  ingest,  calibration,  and  remapping 

2)  navigation 

3)  calibrated  image  retrieval,  magnification, 
geometrical  correction,  and  hard  copy  output 

4)  interactive  analysis  with  ground  truth  data, 
contour  plots,  etc. 

These  four  classes  of  activity  are  necessary  for 
any  satellite  data  analysis  procedure;  similar  systems 
have  been  developed  for  LANDSAT,  SMS,  DMSP 
and  other  geophysical  sensing  satellite  programs.  In 
fact,  as  will  be  noted,  parts  of  the  procedures  to  be 
described  are  adapted  from  other  efforts.  The  pres¬ 


ent  system  is  oriented  towards  polar-orbiting,  rotat¬ 
ing  scanner,  small  acceptance-angle  sensors.  Figure  1 
summarizes  the  overall  procedure. 

Data  Ingest  and  Calibration 

A  qualitative  appreciation  for  passive  remote 
sensing  measurement  techniques  is  necessary  ir 
understand  properly  the  data  ingest  procedure 
design.  Consequently  we  will  give  a  brief  introduc¬ 
tion  to  visible/infrared  remote  sensing  technique 
prior  to  discussion  of  data  ingest  procedures. 

Remote  sensing  systems  make  radiometric  mea¬ 
surements  in  specified  electromagnetic  frequency 
bands  of  energy  emitted  or  reflected  by  the  earth. 
The  absorption  spectrum  for  radiation  passing 
through  the  atmosphere  to  space  from  the  sea  surface 
is  shown  in  Figure  2.  There  are  three  areas  of  this 
spectrum  (“windows”)  in  which  we  have  interest;  the 
visible  band  between  0.3  and  0.7  micrometers,  and 
infrared  bands  between  3.7  and  4.2  micrometers  and 
10.0  and  12.0  micrometers.  These  windows  all  have 
the  characteristic  of  having  relatively  low  absorption 
for  radiation;  thus  they  yield  relatively  clear  views  of 
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Figure  1.  Schematic  view  of  the  satellite  data  processing  pro¬ 
cedure. 
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Figure  2.  Atmospheric  transparency  for  sea  surface  light  as  seen 
from  a  satellite.  Computed  for  a  U.S.  standard  at¬ 
mosphere  utilizing  the  LOWTRAN 4  code  (from  Selby 
et  a!.,  1978). 


the  earth's  surface.  Scanning  radiometers  operating 
in  these  windows  produce  two-dimensional  images  of 
the  sensed  outgoing  energy.  A  point  value  at  the 
spacecraft  for  a  scan-'ine  element  pair  (pixel)  is 
sensed  as  either  a  voltage  or  current  value  propor¬ 
tional  to  the  incident  energy.  This  analog  value  is 
digitized  and  telemetered  to  earth,  or  telemetered  and 
digitized  on  earth.  TIROS-N  is  an  example  of  the 
former*,  and  NOAA-5  of  the  latter'.  Final  digital 
count  values  are  transcribed  into  computer-compati¬ 
ble  format  on  earth.  Possible  error  sources  in  this 
system  occur  during  the  digitization  and  transmission 
steps.  NASA,  NOAA  and  DMSP  spacecraft-data 
streams  contain  known  references  by  which  telemetry 
and  digitizer  errors  can  be  identified  and  corrected. 
Sensor  calibration  is  subject  to  drift  as  well,  due  to 
changes  in  environmental  factors  aboard  the  space¬ 
craft.  This  three-step  process  can  be  described  func¬ 
tionally: 

L  =  L(C(V(A)))  (1) 

where  L  is  a  radiance,  an  energy  measure,  C  a  drift- 
voltage  transfer  function,  V  a  voltage  transfer  func¬ 
tion.  and  A  the  telemetered  signal  level.  In  the  case  of 
infrared  radiometers  the  final  temperature  is  a  non¬ 
linear  function  of  the  radiance  L,  it  is  not  directly 
sensed.  Seawater  can  be  characterized  as  a  grey  body 
emitter  in  these  two  infrared  windows  with  an  emis- 
sivity,  f  of  .98.  To  convert  the  energy  value  to 
temperature,  a  blackbody  transform  is  applied  to  L, 
i.e. 

T  =  T(L/f)  2. 


where  T  is  an  equivalent  blackbody  temperature 
(Planck)  function.  These  transforms  have  been 
documented  by  Koffler*. 

T  is  a  radiometrically  correct  temperature  at  the 
top  of  the  atmosphere.  Two  principal  effects  must  be 
taken  into  account:  variation  in  atmospheric  absorp¬ 
tion  as  a  function  of  scanner  angle  (path  length)  and 
humidity.  Path  length  absorption  is  strongly  con¬ 
trolled  by  the  atmospheric  water  vapor  since  the  ver¬ 
tical  distribution  of  water  vapor  is  not  constant. 
Scientists  have  examined  these  effects  and  shown  that 
the  .sensed  temperature  can  be  up  to  IOC  in  error  for 
very  humid  maritime  atmospheres  at  large  scan 
angles.' "  Two  methods  are  presently  available  in  our 
laboratory  for  calibration.  The  first  is  a  zero-order 
method  based  on  climatological  mean  humidity  pro¬ 
files.  This  was  developed  by  National  Earth  Satellite 
Service  and  is  described  by  Smith'.  A  tempera¬ 
ture  correction  matrix  is  derived  by  regression  as 
function  of  scan  angle  and  sensed  surface  tempera¬ 
ture.  Sensed  surface  temperature  and  total  water 
vapor  concentration  are  found  to  be  highly  cor¬ 
related.  Scenes  with  homogeneous  moisture  fields  are 
readily  analyzed  by  this  technique  which  produces 
semi-quantitative  two-dimensional  surface  tempera¬ 
ture  maps.  Functionally, 

T  =  H(a.T)  (3) 

with  a  the  scanner  angle  id  T  calculated  from  equa¬ 
tion  2.  This  correction,  though  based  on  climatic 
mean  moisture  profiles,  is  ideal  for  quick-look  and 
image  sorting  applications. 
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The  second  method  which  provides  a  belter  cor¬ 
rection  for  atmospheric  absorption  is  based  on  multi- 
spectral  sensing  in  the  water  vapor  bands.  The  new 
AVHRR  radiometers  have  two  or  three  channels  in 
the  infrared:  3.7/r-4.2/i,  IOfr-II;i  and  11/1-12^,  which 
are  channels  3,  4  and  5  respectively.  The  integrated 
atmospheric  transmissivity  over  each  of  these  three 
channels  differs,  and  consequently  algorithms  can  be 
constructed  which  depend  on  the  difference  of  the 
perceived  surface  temperature  in  these  channels.  The 
simplest  algorithm  assumes  that  the  atmosphere  is 
optically  thin  enough  that  the  perceived  surface 
temperature  difference  between  any  of  the  channels 
can  be  parameterized  as  a  simple  function  of  the  bulk 
atmospheric  water  vapor  seen  by  the  satellite  sensor. 
Functionally  a  formula  of  the  following  form  is  ap¬ 
plied; 

Tj  =  ao  +  SLi  Tj  +  aj  (Tj  —  T^)  (3a) 

where  the  Tj’s  are  temperature  retrievals  in  various 
channels  for  a  given  location  and  the  aj’s  give  the 
parameterized  correction.  The  a|’s  can  be  derived 
from  first  principles  theoretically'"  or  can  be  derived 
empirically  from  good  composite  sets  of  surface  and 
satellite  observations."  Both  of  these  methods  give 
better  control  on  the  temperature  accuracy  at  a 
sacrifice  in  processing  throughput. 

Visible  measurements  suffer  from  similar  prob¬ 
lems  although  they  arise  from  different  sources.  An 
analyst  rarely  desires  a  radiance  measurement  at  the 
top  of  the  atmosphere,  but  rather  wishes  to  know  the 
upwelled  radiance  at  the  sea  surface  or  the  mean 
chlorophyll  concentration  or  sediment  load  in  sur¬ 
face  waters.  Whereas  emitted  thermal  radiation  is  an 
intrinsic  property  of  the  medium,  upwelled  visible 
light  arises  from  back-scattered  solar  radiation.  Con¬ 
sequently,  one  must  calculate  the  amounts  of  solar 
radiation  reaching  the  surface  directly  and  scattered 
from  the  atmosphere,  the  resultant  that  is  trans¬ 
mitted  into  the  sea,  that  which  is  back-scattered  by 
the  sea  and  suspended  materials  and  is  transmitted 
through  the  surface  and  back  through  the  atmo¬ 
sphere  to  the  satellite.  The  atmosphere,  ocean  and 
suspended  material  selectively  affect  the  down  and 
upwelled  radiances  as  they  propagate  from  the 
sources  to  the  satellite  sensor.  It  is  a  problem  which 
has  recently  been  straightforwardly  parameterized. 
Atmospheric  correction  algorithms  have  been  pro¬ 
posed.'^  The  method  of  correcting  for  these 
effects'^  and  with  the  clear  water  modifications'*  is 
summarized  in  the  following. 

Radiance  which  reaches  the  sensor  L,(X)  is  com¬ 
posed  of  contributions  from  atmospheric  Rayleigh 


scattering  mX),  aerosol  scattering  L^(X)  and  finally 
the  desired  contribution  L,^(X)  from  the  water.  These 
terms  combine  as: 

L,(X)  =  L,(X)  +  L,(X)  t(X)  L^(X)  (4) 

where  l(X)  is  the  diffuse  transmittance  of  the  atmo¬ 
sphere.  This  formulation  assumes  that  the  sensor  is 
not  viewing  specularly  reflected  sunlight  (glitter)  and 
that  the  fresnel  reflection  from  the  sea  surface  is  in¬ 
cluded  in  the  Rayleigh/aerosol  terms.  The  Rayleigh- 
scattered  light  is  directly  calculable  from  theory. 
Hence,  if  L,^(X)  is  known,  then  L^fX)  is  derivable.  The 
initial  formulation  of  the  correction  process  uses  this 
approach.'’  If  it  is  known  that  the  aerosol  size- 
refractive  index  distribution  is  horizontally  homoge¬ 
neous  over  the  region  of  interest,  then  this  algorithm 
quickly  generalizes. 

This  implies  that  the  aerosol  selectivity: 

S(X„X,)  =  L3(X,)/L3{X,)  (5) 

is  constant  over  the  region  of  interest. 

In  offshore  areas,  such  as  the  Sargasso  Sea  or 
the  Gulf  Stream,  it  is  well  known  that  there  is  very  lit¬ 
tle  terrageneous  material  in  surface  waters;  Figure  3 
shows  that  upwelled  visible  radiance  in  the  red  (670 
nanometer  channel  for  NlMBUS-7)  is  quite  small  as 
a  consequence.  One  thus  can  assert  that  whatever  is 
seen  in  this  part  of  the  spectrum  must  arise  from  the 
atmosphere.  That  is,  we  assume 

i(VLw(\.)  =  0  (6) 

Then  in  a  straightforward  way  one  finds  that 
MX)  =  fl-l(X)]-*(L,(X)  -  L,(X)  - 
S(X,Xo)[L,(Xo)  -  L,{Xo)]l  (7) 

Alternatively,  S(X,  X^)  can  be  computed  with  high  ac¬ 
curacy  if  L„(X)  is  known  for  some  area  in  the  scene. 
Gordon  and  Clark'*  show  that  L^^,(X)  for  Case  1 
waters  can  be  calculated  with  just  a  knowledge  of  the 
solar  zenith  angle  and  that  operationally  these  waters 
are  easily  identified.'’  Thus,  in  this  method,  the  pro¬ 
cedure  is  to  find  a  “clear  water”  area,  compute  L,^(X), 
derive  S(X,  X^),  and  proceed  to  correct  the  image.  A 
case  study  is  presented  later  in  this  report. 

Geometrical  Correction 

Equally  as  important  as  radiometric  correction 
to  the  final  data  product  is  conversion  of  scanner 
elements  to  some  easily  manipulated  aspect  ratio 
(metric).  Some  satellite  systems,  e.g.  DMSP,  ac- 
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complish  a  portion  of  this  in  hardware  by  control  of  pixels  is  defined: 
scanner  or  display  write-rate'’.  NO  A  A  and  most 

other  polar  orbiting  satellites  have  fixed-rate  scan-  F  =  F(s,  e,  <t>,  h)  (8) 

ners  and  are  operated  with  orbits  that  have  inclina-  where  s  is  scan  line,  e  element,  (t>  latitude  and  h 

tions  differing  from  90  degrees*.  Thus,  these  in-  satellite  height.  Figure  4  sketches  the  algorithm.  To 

struments  geometrically  distort  the  pixel  aspect  ratio  first  order,  there  are  two  effects  to  be  corrected:  pixel 

at  larger  scanner  angles.  Other  distortions  due  to  distortion  along  a  scan  line  and  relative  scan  line  mo- 

satellite  vehicle  pitch,  roll  and  yaw  may  also  be  pre-  tion.  The  pixel  distortion  is  taken  to  be  discrete,  i.e., 

sent.  Pitch,  roll  and  yaw  have  much  more  compli-  old  samples  are  mapped  to  the  nearest  whole  new 

cated  interactions  with  the  scan-line  aspect  ratio.  sample.  As  shown  in  Figure  5,  the  vector  contains 

Variations  in  these  parameters  are  deduced  by  study  varying  repetition  factors  “n”,  which  increase  as  a 

of  the  differences  in  high  precision  navigation  function  of  distance  from  the  nadir  point.  A  given 

generated  and  observed  landmark  positions.  old  sample  is  mapped  into  “n”  new  samples.  A  single 

Legeckis  and  Pritchard”  describe  an  algorithm  run-length  encoded  transform  vector  is  constructed 

for  geometrical  correction  of  scan  line  data.  Al-  based  at  the  nadir  point,  since,  to  first  order,  the  cor- 

though  similar  methodology  is  used,  latitude,  time  rection  is  symmetric  about  the  nadir  track.  This 

and  height-dependent  functionals  are  not  based  on  transform  is  .quite  fast  due  to  run-length  encoding  of 

mean  estimates,  but  carry  the  full  dependence  of  each  the  mapping.  The  second  distorting  factor,  earth 

parameter.  Estimates  of  these  parameters  are  avail-  rotation,  to  first  order,  manifests  itself  as  a  line-shift 

able  from  an  ephemeris-driven,  high  precision  navi-  correction.  These  are  either  line  slips  or  advances 

gation  algorithm  that  will  be  described  in  a  latter  sec-  relative  to  the  previous  line  dependent  on  satellite 

tion.  The  geometrical  correction  produces  scan  lines  nadir  track  sense,  scanner  mode  and  hemisphere, 

containing  pixels  with  a  prescribed  aspect  ratio.  Data  sets  that  result  from  transforms  1-8  com- 

Functionally  a  unidirectional  mapping  from  scan-  prise  remapped  scan-line  sets  which  are  stored  on 

element  to  new  constant  aspect  ratio  scan-element  disk  and/or  tape.  These  fixed-format,  calibrated. 


OneM982 


11 


IMAGE  RECTIFICATION 
EQUAL  AREA  MAPPING 

A  »  ?  Ijl  ?  i  I  ^ 

i^3|3l2l2lll2l2l3i^ 


EARTH  ROTATION  CORRECTION 


Figure  4.  Pictorial  represenlaiion  of  scan  tine  and  earth  rotation 
remapping  corrections. 


constant-aspect-ratio  data  sets  constructed  from  the 
polar  orbiter  raw  data,  are  geometrically  and  radio- 
metrically  similar  to  geosynchronous  spin-scan  and 
other  small  acceptance  angle  scanners.  All  subse¬ 
quent  analyses  utilize  this  data-base  store  for  input 
which  simplifies  the  succeeding  levels  of  analysis 
algorithms. 

High  Precision  Navigation 

Orbiting  vehicles  can  easily  be  engineered  with 
accurate  time  bases,  but  it  is  more  difficult  to  main¬ 
tain  precise  vehicle  position  and  attitude  informa¬ 
tion.  Older  sensing  systems,  of  which  NOAA-5  is  a 
prime  example,  have  marginal  attitude  sensors  ( ±  .5 
degree  s'abilization)'.  Newer  vehicles  such  as  TIROS- 
N’,  OMSP-BLOCK  5  D,”  have  accurate  attitude 
( ±  1°  stabilization)  and  temporal  ( ±  50  msec.)  deter¬ 
mination  systems.  The  TIROS-N  time-keeping  ac¬ 
curacy  is  primarily  determined  by  ground-station 
operational  constraints  on  spacecraft  clock  update 


frequency.  Since  a  one-second  error  translates  into 
an  along-track  position  error  of  6.6  km.,  this  is  a 
tight  constraint  on  update  accuracy.  We  observe 
(1979-1980)  variations  as  large  as  two  seconds  (13  km 
along  track  error)  in  clock  accuracy. 

In  all  cases  there  has  been  little  done  to  provide 
useable  onboard  navigation-linked  sensor  process¬ 
ing.  An  ephemeris-driven  algorithm, has  been  con¬ 
structed  to  provide  latitude,  longitude  and  height  in¬ 
formation  for  low-precision  mapping.  A  Keplerian 
model  is  utilized  with  second  and  fourth  order  grav¬ 
ity  terms;  second  order  secular  perturbation  correc¬ 
tions  on  the  period,  mean  anomaly  and  ascending 
node  are  applied  in  addition.  Functionally  the  posi¬ 
tion: 

Px  4h  =  «h  (Ephemeris)  (9) 

is  derived.  The  ephemeris  consists  of  ;he  orbital 
elements:  longitude  of  ascending  node,  angle  of  or¬ 
bital  inclination,  argument  of  perigee,  semi-major 
axis,  eccentricity,  and  universal  time.  X'  is  the 
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longitude,  <t>  latitude,  h  height.  This  procedure  is 
well  known  and  is  lucidly  documented”.  Tests  of  this 
high-precision  navigation  algorithm,  P  demon¬ 
strate  that  satellite  position  can  be  predicted  on  the 
basis  of  daily  mean  orbital  elements  (ephemerides)  to 
±  1  km. 

Precision  Navigation/Mapping 

The  high  precision  of  the  ephemeris-driven 
navigation  cannot  be  fully  utilized  to  establish  precise 
scan,  element  location  due  to  spacecraft  clock  drift 
(c-13  km  induced  error)  and  attitude  control(  ±  4  km 
induced  error  in  roll,  pitch,  and  yaw).  Quantification 
of  these  errors  requires  generation  of  a  precision 
polar-orbit  navigation  model.  Philips  and  Smith” 
describe  a  landmarking  (control  point)  approach  for 
navigation  of  geosynchronous  satellites.  Knowledge 
of  the  scanner  geometry  permits  computation  of  cor¬ 
rections  to  the  ephemerides  by  comparison  of  known 
earth  location  positions  with  model-predicted  posi¬ 
tions.  Daily  update  of  the  model  then  allows  “predic¬ 
tive”  navigation  of  the  spin-scan  radiometer  data. 
This  approach  is  not  directly  applicable  to  the  low 
altitude,  several-pass-per-day  polar-orbiter  situation 
due  to  the  presence  of  non-secular  perturbations. 

A  navigation  method  for  polar  orbiters  has  been 
developed  in  our  laboratory.”  In  this  method,  the 
nadir  track  and  the  scan  lines  are  considered  to  define 
great  circles  on  a  spherical  earth  such  that  the  any 
scan  line  great  circle  is  always  perpendicular  to  the 
nadir  track  great  circle.  (This  is  equivalent  lo  the  ap¬ 
proach  of  Ruff  and  Gruber,”  however,  the  trans¬ 
form  analysis  is  performed  quite  differently.)  The 
procedure  to  establish  navigation  parameters  for  a 
satellite  pass  consists  of  the  following  steps:  First, 
from  the  universal  time  of  the  initial  scan  line,  the 
high  precision  navigation  model  is  used  to  calculate 
the  earth  location  of  the  nadir  pixel.  Second,  land¬ 
marks  (or  continental  outlines)  are  selected  with  an 
emphasis  on  maximizing  their  spatial  spread  with 
respect  to  the  nadir  track.  This  is  done  to  minimize 
the  entry  of  systematic  error  in  the  transform  equa¬ 
tions.  Third,  corrections  from  the  outlines  can  be 
computed  for  the  coordinates  of  the  initial  scan  line, 
the  time  correction  and  mean  pitch,  roll  and  yaw  of 
the  satellite. 

Analytical  transforms  are  generated  by  this  tech¬ 
nique  which  relate  constant-aspect-ratio  pixel  scan¬ 
line  elements  to  latitude,  longitude  and  vice-versa. 
These  transforms  form  a  basis  for  derivation  of 
mapped  projections  of  scanner  data. 

Performance  of  detailed  analyses  of  satellite 
derived  ocean  skin  temperature  also  requires  a 


humanized  interface  for  data  extraction,  enhance¬ 
ment  and  display.  An  intelligent  extraction  algorithm 
founded  on  the  precision  navigation  transforms 
yields  an  easily  useable  method  to  display  an  image 
centered  at  some  earth/scanner  coordinates  with  a 
given  resolution.  This  resultant  data  set  (or  the 
original)  with  an  appropriate  transform  function  can 
then  be  projected  into  various  coordinate  systems 
(transverse  mercator,  polar  stereographic,  etc.) 
which  are  stored  on  disks  and  are  available  for  recall 
at  the  analyst’s  convenience. 

Case  Studies 

In  the  first  part  of  this  report  we  summarized  the 
radiative  transfer,  navigation  and  data  bases  needed 
for  satellite  remote  sensing  of  visible  and  infrared  ra¬ 
diance  from  the  ocean.  Now  we  will  describe  several 
case  studies  which  illustrate  that  application  of  the 
above  methodology  yields  quantitative  and  qualita¬ 
tive  improvement  of  our  understanding  of  the  ocean. 

NOAA-5  Gulf  Stream  Location  Observation 
in  the  South  Atlantic  Bight 

During  early  April,  1977,  Dr.  T.  Lee  (Univ  ersity 
of  Miami/ Rosenstiel  School  of  Marine  and  At¬ 
mospheric  Science)  mapped  the  thermal  structure  of 
the  Gulf  Stream/shelf  water  interface  in  the  Georgia 
Bight.”  A  surface  thermal  field  derived  from  these 
measurements  is  shown  in  Figure  5.  Simultaneously 
the  ocean  skin  temperature  was  observed  from  space 
by  the  NOAA-5  Very  High  Resolution  Infrared 
Radio-meter.  These  data  were  processed  by  the 
afore-mentioned  procedure  and  a  3  km  gridded 
equal-area  product  was  generated  for  analysis  pur¬ 
poses.  This  result  is  illustrated  in  Figure  6.  All  con¬ 
tinental  outlines  and  political  boundaries  have  been 
inserted  from  a  national  data  base  in  terms  of  earth 
location  pairs  via  the  earth-scanner  mar  ing 
transform.  Retrieval  of  the  temperature  values  along 
the  ship  track  is  accomplished  with  the  inverse  mapp¬ 
ing  transform.  The  ship  track  location  and  cor¬ 
responding  temperatures  were  digitized  as  triplets  of 
latitude,  longitude  and  temperature  from  Figure  5 
with  a  graphics  tablet.  Figure  7  shows  the  highlighted 
ship  track  off  the  Georgia  coast.  Magnifying  the  area 
of  interest  about  the  ship  track  (Figure  8),  we  notice 
the  stror  gradients  in  thermal  field  at  the  Gulf 
Stream/shelf  water  interface.  Visual  comparison  of 
Figures  5  and  7  illustrates  the  fact  that  the  satellite 
data  with  its  inherently  higher  spatial  definition  can 
yield  a  good  sarface  thermal  map  in  situations  with 
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Figure  6.  NOAA  S  VHRR  equal-area  processed  thermal  map  for 
10  April  1977. 


low  to  moderate  atmospheric  water  vapor  (<  1 
gm/cm’)-  [A  mean  airmass  correction  of  1.25  C  is 
subtracted  from  the  satellite  observations.]  The 
temperature  along  the  ship  track  comprises  a  one 
dimensional  series  as  a  function  of  cumulative  along- 
track  distance.  Frontal  definition  and  temperature 
agreement  between  both  methods  is  quite  good. 
Quantitatively  we  find  the  standard  deviations  of  the 
ship  and  satellite  track-line  series  to  be  2. 1C  and 
2.2C,  respectively,  while  the  rms  difference  of  the 
two  series  is  1 .2C.  These  values  are  calculated  for  a 
set  of  several  thousand  observation  pairs.  The  near 
equivalence  of  the  ship  and  satellite  standard  devia¬ 
tions  implies  there  is  little  second  order-moment  bias, 
or  conversely  the  satellite  shows  the  same  biases  as 
the  ship.  A  mean  square  difference  of  1 .2C  is  the  ex¬ 
pected  difference  at  ±  1  standard  deviation.  This  is 
comparable  to  estimates  NOAA/NESS  has  made  for 
other  temperature  products.”  We  believe  the  main 
contributor  to  this  difference  is  the  inherently  dif¬ 


ferent  sampling  scales  of  the  two  methods.  The  mean 
square  and  mean  differences  are  extraordinary  in 
that  the  comparison  values  are  representative  of  total 
error,  that  is  errors  arising  from  navigation  inac¬ 
curacies,  time  and  space  sampling  mismatch,  water 
vapor  contamination,  radiometer  calibration,  etc. 
This  verification  area  was  chosen  for  two  reasons: 
high  quality  ship  and  satellite  data  were  available 
nearly  concurrently  and  strong  thermal  features  with 
small  scales  were  present  which  would  magnify 
navigation,  mapping,  and  calibration  errors.  Thus 
we  believe  that  this  comparison  is  indicative  of  the 
procedure’s  utility  in  relatively  clear,  homogeneous 
atmospheric  water  vapor  situations.  A  spatially  in¬ 
homogeneous  moisture  field  would  render  this 
analysis  less  valid. 
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Figure  7.  Magnification  of  Figure  6  about  the  ship  track. 


NIMBUS-7  CZCS  Observations  of  Color 
Variability  in  the  Middle  Atlantic  Bight 

A  recent  work  by  Gordon  et  aP*  utilizes  CZCS 
(Coastal  Zone  Color  Scanner)  observations  to  map 
the  surface  chlorophyll  structure  of  a  warm-core  ring 
off  the  Northeast  coast  of  the  United  States.  Warm- 
core  rings  occur  when  extended  meanders  of  the  Gulf 
Stream  pinch  off  on  the  cold-wall  (north  side)  and 
form  a  closed  anticyclonic  rotating  vortex  of  Gulf 
Stream  and  Sargasso  Sea  waters.  This  vortex  of  Gulf 
Stream  and  Sargasso  Sea  waters  is  imbedded  in  the 


surrounding  slope  or  shelf  waters  and  is  thought  to 
isolate  an  ecosystem  at  some  point  in  its  life  cycle. 
Rings  have  horizontal  dimensions  of  l(X)-200km  and 
vertical  extents  equal  to  or  greater  than  1000m.-'' 
We  show  in  this  study  that  mean  chlorophyll  concen¬ 
tration  in  surface  waters  can  be  observed  to  35  “/o 
utilizing  satellite  method  alone. 

Figure  9  is  an  overview  of  the  shelf  and  slope 
water  chlorophyll  concentration  off  the  northeast 
United  States  coast  for  15  June  1981 .  The  readers  at¬ 
tention  is  directed  to  several  prominent  features:  the 
low  chlorophyll  values  (lighter  grays)  in  a  warm  core 
ring  located  at  39°  N,  69°  W,  (catalogued  as  79-B  by 
Fitzgerald  and  Chamberlin”)  the  low  and  high  chlo¬ 
rophyll  streamers  wrapping  into  the  ring,  their 
length,  and  the  overall  decrease  of  chlorophyll  as  one 
travels  north  of  the  Gulf  Stream  and  away  from  the 
coast. 


INTERCOMPARISON  OF  SHIP  VS.  SATELLITE  SST  FOR  04/10/77 

Figure  8.  Satellite  versus  ship  derived  temperature  intercomparison. 
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The  fact  that  this  warm  core  ring  has  such  a 
strong  chlorophyll  signature  implies  that  it  can  be 
tracked  as  either  a  color  or  an  infrared  anomaly. 
Qualitatively  the  color  structure  raises  several  in¬ 
teresting  theoretical  questions:  How  is  it  that  the  sur¬ 
face  chlorophyll  streamers  cross  the  ring-front  and 
spiral  in  the  central  regions  of  the  ring?  What  is  the 
mechanism  underlying  propagation  of  ring  effects 
(through  the  streamers)  to  distances  of  one  to  two 
ring-diameters  away?  Are  there  differences  in  the 
thermal  and  color  signatures  of  a  warm-core  ring? 
Does  a  regional  water  mass  analysis  support  the 
large-scale  chlorophyll  surface  pattern  we  observe? 
Comprehensive  delineation  of  the  first  two  effects 
awaits  1982  warm  core  ring  cruise  work,  the  latter 
question  can  be  addressed  qualitatively  by  looking  at 
the  slope  water  mass  distributions  in  the  area. 
McLellan’s  analysis"  of  the  water  mass  areal  distri¬ 
butions  in  this  area  agrees  with  the  present  observa¬ 
tions  in  that  he  finds  shelf-type  waters  in  the  Bight 


the  different  sensing  methods,  however  delineation 
of  cross  boundary  flows  is  quite  a  bit  more  difficult 
in  the  infrared.  The  local  minimum  in  central  waters 
seen  in  the  chlorophyll  is  mirrored  in  the  thermal  im¬ 
agery  as  a  temperature  maximum.  One  is  left  with  the 
impression  that  although  there  is  coherency  between 
the  two  methods,  color  is  a  more  sensitive  indicator 
for  study  of  cross-boundary  surface  flows. 

Fine-scale  features  may  be  atmospheric  arti¬ 
facts.  One  must  know  which  featuris  are  atmo¬ 
spheric  and/or  oceanic  in  origin.  A  cru  se  in  support 
of  the  NIMBUS  CZCS  validation  made  a  section 
through  this  ring  at  0000  EDT,  10  June  1979.  NIM¬ 
BUS  orbits  3157  and  3171  bracket  this  ship  coverage 
occurring  approximately  12  hours  before  and  after, 
respectively,  ship  passage  through  ring  center. 

Pigment  concentration  was  measured  along  the 
ship  track  by  using  a  Turner  Designs  flow-through 
fluorometer  (Model  10)  to  monitor  the  in-vivo  fluo¬ 
rescence  of  water  drawn  continuously  from  about 


and  up  into  the  Gulf  of  Maine  and  cooler,  lower  ox¬ 
ygen  waters  of  Labrador  origin  offshore.  This  jux¬ 
taposition  of  water  masses  yields  the  strong  horizon¬ 
tal  chlorophyll  contrasts  seen  in  the  present  imagery. 

Figure  10  is  an  enlarged  view  of  Figure  9 
centered  about  the  ring.  Streamer  delineation  and 
local  chlorophyll  minimum  in  the  ring  interior  are 
seen.  Figure  11  is  an  enlarged  view  of  the  ring  as  seen 
several  hours  earlier  in  TIROS-N  AVHRR  infrared 
imagery.  Ring  frontal  outlines  are  consistent  between 


Figure  9.  NIMBUS-7  CZCS  derived  chlorophyll  image  for  15 
June  1979,  orhil  3226.  Lighter  shades  denote  lower 
values,  darker  higher  (from  Gordon  et  ai,  1982). 

1.5m  below  the  surface.  Periodically,  discrete 
samples  were  drawn  from  the  same  source  and  ex¬ 
tractions  performed  to  determine  fluorometrically 
the  concentrations  of  chlorophyll  and  the  associated 
phaeopigments,  using  the  technique  described  by 
Yentsch  and  Menzel”  with  the  modifications  of 
Holm-Hansen  et  aP*.  These  discrete  determinations 
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were  then  used  to  relate  the  continuous  in-vivo  fluo¬ 
rescence  measurements  to  the  pigment  concentration. 
CZCS-derived  pigment  concentration  (heavy  line)  is 
compared  with  these  continuous  determinations  in 
Figure  12.  This  track  is  from  south  to  north,  starting 
in  the  Gulf  Stream,  passing  through  the  slope  water, 
the  warm-core  ring  and  finally  ending  on  the  shelf. 
These  distinct  water  types  are  apparent  in  Figure  12. 
The  cold  watl  of  the  Gulf  Stream  is  at  about  220km 
on  the  axis,  and  the  ring  extends  from  300  to  350km. 
Satellite  determined  pigment  concentration  is  split 
between  orbit  3157  from  320  to  440km  and  orbit  3171 
from  420  to  600km.  Gaps  in  the  satellite-derived 
values  are  due  to  clouds.  Agreement  between  the  two 
methods  is  apparent. 

Examination  of  warm-core  rings  is  an  ideal  ap¬ 
plication  for  remote  sensing  since  they  have  signa¬ 
tures  in  the  color  and  infrared.  Technically  they  are 
easy  to  study  as  well  since  the  Gulf  Stream  offers  a 
dear-water  area  for  CZCS  calibration  and  slowly 
varying  thermal  signal  for  infrared  observations. 


TlROS-N/NOAA-6  Observations  of 
Somali  Current  Response  to 
the  Summer  Monsoon 

Our  purpose  in  this  section  is  to  give  a  summary 
of  surface  temperature  variability  near  the  East 
African  Coast  as  observed  by  satellite  and  ship  dur¬ 
ing  the  mature  phase  of  the  summer  monsoons  of 
1976-1979.  Sea  surface  temperature  (SST)  pattern 
changes  have  been  shown  to  be  strongly  correlated 
with  changes  in  the  surface-current  field”  as  such 
they  yield  a  sequence  of  observations  of  the  East 
African  current-system  variability  during  the  summer 
monsoon. 

Monsoonal  forcing  has  two  major  large-scale  ef¬ 
fects  on  the  northwestern  Indian  Ocean:  Somali  Cur¬ 
rent  spinup  and  maintenance,  and  cooling  of  the  up¬ 
per  mixed  layer.  Delineation  of  these  effects  on  the 
ocean  is  complicated  by  the  presence  of  large  scale 
eddies  near  the  ocean  boundary  (which  maybe  the 
Somali  Current)  and  in  the  ocean  interior”,  intensive 
upwelling  near  the  boundaries”,  and  other  large  scale 
disturbances. 

Historically  physical  oceanographers  have  con- 
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Figure  12.  Comparison  of  ship  and  satellile  derived  surface  pig¬ 
ment  concentrations  for  NIMBUS-7  orbits  3157  and 
3171,  and  R/y  Athena  II  in-situ  observations.  Ship 
observations  are  denoted  by  the  light  line,  satellile 
observations  by  the  heavier  solid  and  dashed  lines,  or¬ 
bits  3157  and  3171  respectively  (from  Gordon  el  al., 
1982). 


sidered  the  Somali  Current  to  be  a  seasonally  revers¬ 
ing  western  boundary  current.  Source  waters  from 
south  of  the  equator  flow  north  along  the  East 
African  Coast,  separating  and  flowing  into  the  inte¬ 
rior  as  part  of  the  “Great  Whirl”  between  10"  and 
“  the  Great  Whirl  being  a  large  anti- 
cyclonic  eddy  which  is  observed  off  the  coast  of 
Somalia  during  the  summer  monsoon.  It  has  dimen¬ 
sions  of  500  km  (parallel  to  coast)  by  400  km  (perpen¬ 
dicular  to  the  coast. “  Bruce  and  Brown  find  several 
instances  (1970,  1976,  1979)  when  the  Somali  Current 
is  observed  to  separate  from  the  coast  south  of  the 
Great  Whirl  between  latitudes  2°  N-5°  N.“  ”  It  was 
assumed  that  there  were  “single  gyre  years”  and 
“two  gyre  years.”  In  1979  Brown  found  both  cases  at 
different  stages  of  the  summer  monsoon.” 

Figure  13  illustrates  the  ship  derived  surface 
thermal  fields  as  observed  during  the  mature  phase  of 
the  summer  monsoon  in  1979”.  In  the  surface- 
current  field  at  this  time  offshore  flow  was  found  at 
2  "  N,  together  with  an  anti-cyclonic  gyre  between  5  ° 
N  and  10°  N  with  offshore  flow  from  10°  N  to  12° 
N.“  Surface  current  and  thermal  fields  are  strongly 
coupled  as  might  be  expected  from  earlier  observa¬ 
tions”.  Separation  regions  are  expressed  in  the  SST 


as  cold  wedge  shaped  areas. 

Figure  14  presents  a  time  series  of  frontal  posi¬ 
tions  as  observed  during  the  period  June  to  October, 
1979,  by  satellite.  The  southern  frontal  region  starts 
to  move  poleward,  slows  and  then  quickly  translates 
to  a  position  ~  1  °  south  of  the  original  northerly 
separation  location.  During  the  same  time  interval 
the  northerly  region  (dashed  area)  moves  slowly 
poleward  to  the  region  of  Socotra  and  thence  more 
quickly  southward,  coalescing  with  the  southern 
region  south  of  its  earlier  position  near  Ras  Hafun. 
Ship  observations  made  during  the  latter  phases  of 
this  process  near  Ras  Hafun  find  fresher,  cooler 
water  in  the  upper  - 100  m.  overlaying  warmer, 
more  saline  waters  in  the  northwest  quadrant  of  the 
Great  Whirl.  That  is,  surface  waters  are  characteris¬ 
tic  of  Somali  Current  origin  while  deeper  waters  are 
compatible  with  Arabian  Sea  sources  (John  Bruce, 
private  communication.  Woods  Hole  Oceanographic 
Institution). 

Analysis  of  the  1979  results  stimulated  interest  in 
reprocessing  of  earlier  satellite  data  archived  during 
1976-1978.  Figure  15  illustrates  observed  frontal  time 
series  for  1976.  Southern  frontal  translation  occurs 
later  and  is  much  slower  than  the  1979  case.  The 
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Figure  13.  Time  sequence  of  ship  observed  sea  surface  temper¬ 
ature  fields  for  oceanic  response  to  the  southwest  mon¬ 
soon,  1979.  la:  16  April-10  May,  lb:  II  May-22  May, 
Ic:  23  Muy-12  June,  Id:  22  June-31  July,  le:  18 
Augusl-3  September,  (from  Brown  et  al,  1980). 


multi-phase  nature  of  the  transition  is  evident,  al¬ 
though  the  absolute  changes  in  translation  rate  are 
much  less. 

Satellite  imaging  of  the  area  in  1977  was  inade¬ 
quate  for  resolution  of  weaker  thermal  differences 
due  to  spacecraft  malfunctions  and  related  problems. 
Imagery  on  hand  neither  supports  nor  denies  the 
presence  of  the  dual  separation  region  in  this  year. 
Research  ship  observations  of  surface  temperature 
are  available  for  1977  from  the  R/V  Shokals’kii 
which  was  deploying  INDEX  moorings  off  Somalia. 
A  temperature  map  for  early  June  is  presented  in 
Figure  16  (Donald  Johnson,  private  communication. 
Old  Dominion  University,  Norfolk).  Temperatures 


are  not  as  low  as  in  the  other  three  years,  but  a  wedge 
is  apparent  at  6  N. 

Satellite-derived  frontal  positions  during  sum¬ 
mer  1978  are  illustrated  in  Figure  17.  Fronts  occur  at 
the  same  times  as  in  1979,  however,  the  northern 
frontal  region  passes  through  the  Socotra  topo¬ 
graphic  arc  and  eventually  dissipates  in  the  western 
Arabian  Sea.  Southern  wedge  translation  rates  are  in¬ 
termediate  between  those  observed  in  1976  and  1979. 

Satellite-derived  frontal  patterns  suggest  that 
dual  separation  regions  (double  gyre  configuration) 
occur  in  three  of  the  four  years  studied.  This  inter¬ 
pretation  is  at  variance  with  classical  views  of  Somali 
Current  circulation. 
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Figure  16.  Surface  temperature  map  derived  from  measurements 
made  onboard  the  R/V  Shokais'kii  during  May-June, 
1977  off  Somalia,  (from  Evans  and  Brown.  1981, 
courtesy  Dr.  D.  Johnson). 


Figure  14.  Time  .sequence  of  satellite  observed  frontal  locations 
for  northern  and  southern  frontal  wedges  for  year 
1979.  The  northern  wedge  front  is  hatched  (from 
Evans  and  Brown,  1981).** 
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Figure  15.  Time  sequence  of  sateldte  observed  frontal  locations 
for  northern  and  southern  frontal  wedge'i  for  year 
1976.  The  northern  wedge  front  i.s  hatched  (from 
Evans  and  Brown.  1981). 


These  observations  show  coalescence  of  the 
Somali  Current  system  is  abrupt  and  can  be  thought 
to  consist  of  four  phases.  During  phase  I  the 
southern  separation  region  moves  north  from  its 
onset  location  during  the  mature  phase  of  the  sum¬ 
mer  monsoon  to  -6°  N.  Simultaneously  the  north¬ 
ern  frontal  region  translates  poleward  to  the  region 
of  Socotra.  In  phase  II  the  system  is  in  quasi-steady 
state  for  a  period  of  several  days  to  a  week.  The  third 
phase  begins  with  reversal  of  northern  frontal  propa¬ 
gation  in  1976  and  1979;  in  1978  phase  III  is  initiated 
when  the  northern  frontal  zone  breaks  through  the 
topographic  arc  connecting  Socotra  and  East  Africa. 
The  southern  region  moves  more  quickly  poleward 
than  in  phase  I,  while  the  northern  zone  translates 
equatorward  (1976,  1979);  in  1978  the  northern  zone 
translates  into  the  Gulf  of  Aden  and  dissipates. 
Finally  coalescence  occurs  during  phase  IV  in  1976 
and  1979.  This  phase  lasts  for  periods  of  days  to 
weeks. 

From  the  data  presented  for  the  period  1976- 
1979  we  assert  that  this  “anomalous”  separation  is 
confirmed  in  3  of  the  4  years.  More  interestingly 
these  satellite  data  give  us  a  first  look  at  the  temporal 
evolution  of  the  Somali  Current  from  the  anomalous 
separation  case  to  the  more  historically  expected  cir¬ 
culation.  It  may  well  be  that  the  single  or  dual  gyre 
nature  historically  observed  is  a  function  of  the  time 
when  observations  were  made  rather  than  being  in¬ 
dicative  of  interannual  differences  in  oceanic 
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Figure  17.  Time  sequence  of  satellite  observed  frontal  locations 
for  northern  and  southern  frontal  wedges  for  year 
1978.  The  northern  wedge  front  is  hatched  (frorn 
Evans  and  Brown,  1981). 


response. 

Spinup  model  results  suggest  that,  at  least  in 
some  years,  the  Somali  Current  can  separate  from 
the  coast  south  of  the  Great  Whirl  at  5°  N.  This 
variability  may  be  due  to  the  strength  of  the  mon¬ 
soon.  Cox  has  shown  the  generation  and  stall  char¬ 
acteristics  of  the  eddy  circulations  are  strongly 
dependent  on  the  the  wind  stress. “  If  one  assumes 
that  there  is  a  transition  from  a  stalled-eddy  state  to  a 
poleward-translating  state  with  diminution  of  the 
wind  stress,  then  the  Cox  (1979)  model  may  also  be 
relevant  to  the  phase  I  observations. 

Conclusions 


of  SST  response  to  the  summer  monsoon?  This  is  not 
to  say  that  quantitative  results  are  unnecessary;  they 
directly  give  us  an  appreciation  for  atmospheric  con¬ 
tamination  of  retrieval  and  enable  us  to  judge  which 
signals  are  oceanic  in  origin.  It  is  important  to 
understand  that  these  types  of  passive  measurements 
permit  generalization  of  accurate  in-situ  process 
observations  to  larger  scales,  but  do  not  at  this  time 
have  the  accuracy  of  in-situ  observations.  Satellite 
oceanography  currently  is  a  highly  descriptive  en¬ 
deavor,  only  in  limited  cases  will  it  be  possible  to  ob¬ 
tain  very  quantitative  results.  However  even  in  its 
presently  crude  quantitative  form  it  does  and  will 
enable  workers  to  test  hypotheses  of  ocean  circula¬ 
tion  and  transient  response  to  forcing. 
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Introduction 

The  great  difficulty  of  acquiring  requisite 
observational  data  plagues  all  efforts  at  ocean 
analysis  and  forecasting.  In  recent  decades  ocean¬ 
ographers  and  acousticians  have  acquired  measure¬ 
ments  that  indicate  kinetic  energy  at  mesoscales  (< 
24  hours  and  tens  to  hundreds  of  kilometers)  often 
exceeds  that  of  the  synoptic  mean  motion  by  as  much 
as  two  orders  of  magnitude.'  The  importance  of 
mesoscale  processes  and  the  means  to  observe  and 
model  them  now  commands  the  attention  of  much 
modern  ocean  research.  Our  understanding  of  meso¬ 
scale  processes  in  coastal  or  oceanic  boundary 
regions  has  at  present  inadequate  observational  and 
theoretical  bases.  Oceanic  processes  at  those  scales 
generally  defy  prediction  or  even  detection.  Satellite 
scanning  radiometers  have  the  potential  to  supply  the 
needed  synoptic  data  but  for  two  serious  problems. 
First,  the  atmosphere,  even  without  clouds,  can 
obliterate  the  desired  oceanographic  information. 
Secondly,  the  radiometric  temperatures  may  repre¬ 
sent  only  the  skin  temperature  of  the  sea.  We  have 
found,  and  describe  here,  those  times  and  places 
where  these  problems  are  minimized  and  a 
“window”  on  the  mesoscale  structure  of  the  upper 
ocean  occurs. 

Over  the  continents  in  winter,  the  vast  high- 
pressure  systems  of  cold,  dry,  and  dense  polar  air  are 
generated  largely  by  radiative  ooling.  The  mobility 
of  weather  systems  in  mid-latitudes  results  in  periodic 
outbreaks  of  this  air  across  coastal  regions  into  the 
adjacent  ocean  basins.  The  “window”  results  from 
two  aspects  of  this  process.  First  is  the  development 
of  clear,  dry,  relatively  uniform  air  over  the  sea, 
through  which  best  quality  radiometric  measure¬ 
ments  can  be  made.  Second  is  the  deeply  mixed  con¬ 
ditions  caused  by  strong  winds  and  negative  oceanic 
heat  fluxes,  which  make  surface  measurements  repre¬ 
sentative  of  much  of  the  upper  water  column.  Under 


these  conditions  infrared  and  visible  sensors  of  satel¬ 
lites  can  detect  the  mesoscale  patterns  of  surface  tem¬ 
perature  and  turbidity  particularly  well.  These  out¬ 
breaks  occur  in  3-6-day  cycles,  bringing  one  or  more 
days  of  ideal  imaging  conditions  at  a  time  and  yield¬ 
ing  imagery  of  excellent  quality. 

Atmospheric  Effects 

Implicit  in  the  above  discussion  is  the  fact  that 
polar  continental  air  is  superior  to  maritime  air  for 
satellite  observations  and  measurements.  Satellite 
scanning  radiometers  collect  the  planetary  radiances 
that  reach  altitudes  of  “inner  space”  (850  km  for 
NOAA  satellites).  The  radiances  in  the  thermal  in¬ 
frared  part  of  the  spectrum  are  a  combination  of 
radiant  energy  emitted  from  the  earth,  cloud,  or  sea 
surface  that  has  been  altered  more  or  less  by  the 
intervening  atmosphere.  This  alteration  by  a  cloud- 
free  atmosphere  is  a  complex  function  of  molecular 
absorption  and  reemission  processes.  In  the  simplest 
terms,  a  variable  fraction  of  the  surface-emitted  in¬ 
frared  radiation  is  absorbed  and  the  energy  is  re¬ 
emitted  at  a  different  temperature,  mainly  by  the 
triatomic  molecules  of  the  atmosphere  (HjO,  O3,  and 
CO2).  Of  these,  water  vapor  is  the  most  effective  and 
variable  in  distribution.  Its  effect  varies  along  the 
radiation  path  as  a  function  of  its  concentration  and 
temperature.  Under  warm  and  moist  conditions  the 
sea-surface  temperature  patterns  may  be  completely 
obliterated  even  without  condensation  and  cloud  for¬ 
mation.  The  entire  absorption  and  reemission  effect 
is,  however,  minimal  in  winter  when  cold,  dry  air 
overlies  the  lands  and  seas  of  mid-  and  high-latitude 
regions. 

The  whole  problem  of  attenuation  of  infrared 
measurements  has  been  approached  by  most  investi- 
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gators  from  the  standpoint  of  correcting  satellite- 
derived  effective  radiation  temperatures  to  closely 
approach  true  surface  temperatures  measured  dir¬ 
ectly  by  ships.  Through  all  types  of  atmospheres,  in¬ 
vestigators  have  found  satellite  temperatures  lower 
by  1-10°C.  A  uniform  cooling  of  the  field  of  sea- 
surface  temperatures  is  not  so  serious  a  problem  for 
sea-going  users.’  The  satellite  data  can  be  calibrated 
with  surface  measurements  transmitted  from  ships  or 
buoys  to  make  full  use  of  the  finely  resolved  horizon¬ 
tal  patterns  of  temperature.  Cooling  of  radiation 
temperatures  is,  however,  not  the  whole  attenuation 
problem.  Our  studies  of  the  satellite-  and  aircraft- 
based  surface  temperature  measurements  over 
coastal  regions  indicate  that  the  atmospheric  effects 
can  be  usefully  separated  into  two  problems; 

(1)  Atmospheric  attenuation — the  lowering  of 
radiation  temperatures  from  the  surface- 
measured  values  due  to  the  imperfect  trans¬ 
parency  of  the  atmosphere. 

(2)  Suppression  of  temperature  differences — the 
reduction  of  measured  surface  temperature 
differences  between  warm  and  colder  parts 
of  the  scene  caused  by  the  wavelength  varia¬ 
tion  of  radiative  transfer  through  the  at¬ 
mosphere. 

This  second  problem,  from  a  practical  standpoint,  is 
more  serious,  since  it  means  that  atmospheric  effects 
could  obliterate  oceanic  fronts  and  mesoscale  fea¬ 
tures  of  tactical  interest,  to  mislead  the  user. 

A  number  of  studies  of  atmospheric  effects  on 
infrared  radiation  measurements  of  sea-surface 
temperature  have  been  conducted  at  aircraft  alti¬ 
tudes.  In  June  1971  off  the  Massachusetts  coast  an 
aircraft  radiometer  was  unable  to  detect  temperature 
changes  of  2°C  in  the  sea  surface  with  a  moist  layer 
of  about  300  m  depth  present  over  the  sea.’  Another 
well-documented  case  was  reported  by  scientists  of 
the  Research  Division,  Japan  Research  Committee  of 
Environmental  Remote  Sensing  from  aircraft 
multispectral  scanner  flights  over  Ofunato  Bay,  Hon¬ 
shu,  Japan.  In  the  morning  on  25  July  1978  the  in¬ 
frared  channel  showed  excellent-quality  imagery, 
with  absolute  temperatures  slightly  cooler  than 
ground  truth,  as  expected.  That  afternoon,  thermal 
contrast  was  almost  nil,  and  aircraft  measurements 
of  surface  temperatures  were  warmer  than  surface 
truth  measurements.  A  review  of  the  experimental 
results  revealed  that  the  morning  flights  had  been 
conducted  after  the  nighttime  land  breezes  had 
brought  drier,  cooler  air  out  over  the  coastal  bay.  In 
the  afternoon  the  sea  breeze  had  brought  warm, 
moist  air  in  from  the  Pacific,  which  obliterated  the 


gradients  in  the  same  manner  as  we  have  seen  in 
satellite  data.  Radio-sonde  measurements  available 
simultaneously  with  the  flights  show  that  the  2  km- 
thick  atmospheric  column  beneath  the  aircraft  of  the 
morning  flight  was  an  average  of  2  °C  cooler  and  had 
I  g  cm  -  less  precipitable  water  than  the  afternoon 
air  mass. 

Atmospheric  suppression  of  surface  temperature 
differences  is  well  illustrated  with  the  NOAA-4  in¬ 
frared  image  (10.5-12.5  fim  band)  acquired  at  OlOOZ 
on  4  April  1976  over  the  Florida  Peninsula  (Figure  1 ). 
The  temperature  gradients  associated  with  the  Gulf 
Stream  are  seen  to  disappear  near  the  bottom  of  the 
image.  The  Florida  Peninsula  at  this  time  was  partly 
overrun  by  an  atmospheric  cold  front  moving  rapidly 
south.  The  air  mass  boundary  beiween  the  continen¬ 
tal  air  and  the  tropical  maritime  air,  visible  just  south 
of  Lake  Okeechobee,  was  unusually  cloud  free.  Sur¬ 
face  observations  indicate  zero  cloud  cover  for  all 
stations,  except  for  two-tenths  registered  at  West 
Palm  Beach  on  the  east  coast,  which  was  directly  in 
line  with  the  clouds  seen  seaward  on  the  image.  This 
cloud  mass  is  the  landward  extension  of  the  cloud 
shield  of  the  cold  front  visible  over  the  Atlantic.  Air 
temperatures  are  20-21  °C  with  dew  point  readings  of 
2-4 “C  just  north  of  the  front,  and  26 °C  with  dew 
point  readings  of  18-20°C  to  the  south.  There  is  a 
35*7o  humidity  drop  and  a  10  g  m '  moisture  content 
difference  in  the  air  between  the  surface  stations 
across  the  front  (Vero  and  West  Palm  Beach 
stations).  Analysis  of  the  Miami  radiosonde  data 
before  and  after  the  frontal  passages  indicates  that 
the  warm,  moist  atmospheric  conditions,  with  2.6-3 
gcm“^  of  total  precipitable  water  (TPW),  were 
replaced  by  colder,  dry  air  with  1.1-1.5  g  cm"-  TPW. 
Nearly  all  surface  temperature  differences  at  the  edge 
of  the  Gulf  Stream  and  Florida  Peninsula  have  been 
suppressed  in  the  radiation  temperature  field  by  these 
atmospheric  conditions. 

Maul  and  Sidran*  calculated  that  an  8°C 
temperature  difference  at  the  sea  surface,  when 
viewed  vertically  through  the  30 °N  July  Standard  At¬ 
mosphere,  would  be  recorded  as  less  than  3°C  by 
satellite  and  further  reduced  to  1  °C  when  viewed 
through  the  atmosphere  at  an  angle  of  60°  off  the 
vertical.  Using  the  Maul  and  Sidran  model,  the  radi¬ 
ative  transfer  process  through  the  two  air  masses  over 
Florida  at  the  time  of  the  image  (Figure  1)  has  been 
simulated,  and  the  results  are  shown  in  Figure  2.  The 
solid  line  is  the  distribution  of  atmospheric  water 
vapor  in  terms  of  absolute  humidity  versus  altitude. 
The  dashed  lines  are  plots  of  the  calculated  radiation 
temperatures,  as  they  would  be  measured  by  a  down¬ 
ward-looking  radiometer,  rising  through  the  atmo- 
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sphere  for  surface  radiation  temperatures  of  20  and 
22  °C.  Starting  with  the  2°  temperature  difference  at 
the  sea  surface  (20  and  22  °C),  both  atmospheres 
reduced  the  temperatures  and  the  temperature  dif¬ 
ferences  calculated  for  the  top  of  the  atmosphere. 
The  warm,  humid  air  reduced  the  temperature  dif¬ 
ference  from  2  to  0.8  °C,  whereas  the  coolei ,  dry  air 
reduced  the  temperature  difference  from  2  to  1.4°C. 
Examination  of  the  moist  air  plot  shows  that  the  sup¬ 
pression  of  the  temperature  difference  has  occurred 
in  the  lower  moist  boundary  layer  of  the  atmosphere 
(lower  3  km).  The  calculated  radiation  temperature 
has  been  increased  in  the  lower  3  km  prior  to  the 
decrease  in  the  upper  atmosphere.  This  is  a  result  of 
the  infrared  radiance  contribution  made  to  the  over¬ 
all  upward  radiation  by  the  warm,  humid  maritime 
air.  The  suppression  of  temperature  differences  is  the 
cumulative  effect  of  two  processes.  The  Lambert- 
Beer  Law  states  that  the  decrease  in  spectral  intensity 
of  radiation  by  absorption  with  distance  is  a  function 
of  the  intensity  of  the  initial  radiation  as  well  as  of 
the  concentration  of  the  absorbers.  The  second  proc¬ 
ess  is  the  emission  of  infrared  energy  by  the  atmo¬ 
sphere  itself.  The  result  of  both  processes  is  to  reduce 
the  effective  brightness  temperature  of  a  warmer 

Radiation  temperature  (K) 

290  292  294  296  298 


Absolute  humidity  |g/m3) 


body  more  than  that  of  a  cooler  body. 

Our  field  data  suggest  that  atmospheric  attenua¬ 
tion  occurs  under  all  conditions  but  that  the  suppres¬ 
sion  of  temperature  differences  is  an  air-mass- 
dependent  phenomenon.  We  have  attempted  to 
quantify  this  effect  by  obtaining  surface  temperature 
measurements  of  two  distinct  bodies  of  water  simul¬ 
taneously  with  an  offshore  atmospheric  sounding 
and  a  satellite  overpass.  We  accomplished  this  five 
times  with  the  NOAA-5  satellite  in  February  1977  in 
the  coastal  waters  of  the  north-eastern  Gulf  of  Mex¬ 
ico.’  In  the  cases  of  higher  levels  of  TPW,  the  ther¬ 
mal  contrast  between  the  warm  and  cold  waters  was 
suppressed,  and  the  suppression  appears  to  increase 
with  increasing  TPW.  In  the  cases  of  cold,  dry  air 
over  the  sea,  the  atmospheric  effect  appears  to  be 
primarily  one  of  uniform  lowering  of  surface 
temperatures,  without  loss  of  gradients  in  the  data 
field.  The  few  observations  available  with  the  older 
lower  resolution  satellite  radiometer  provide  only  a 
preliminary  look  at  the  “water-vapor  versus 
infrared-imagery  contrast”  relationship  which  re¬ 
quires  further  study.  A  number  of  atmospheric 
radiative  transfer  models  are  available’  *,  but  we 
have  found  that  they  seriously  overestimate  the  ai- 
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Figure  2.  Atmospheric  humidity  profiles  and  simulated  radiative  transfer  of  in  frared  energy  for  20  and 
22°C  water  surface  through  the  two  air  masses.  S  and  4  April  1975. 
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mospheric  attenuation.  They  do  somewhat  better  in 
estimating  temperature  difference  suppression. 

Some  Mesoscale  Structures  of 
the  Upper  Ocean 

Examples  of  the  best  quality  infrared  imagery 
available  from  the  new  generation  of  NOAA 
satellites  are  drawn  from  the  study  of  the  Korea 
Strait  and  adjacent  seas  sponsored  by  the  ONR 
Coastal  Sciences  Program.  Perhaps  the  most  severe 
outbreaks  of  polar  continental  air  in  the  world  occur 
over  this  oceanic  region  off  northeast  Asia.  A  great 
variety  of  marine  environments  occur  here,  ranging 
from  the  Kuroshio-dominated  edge  of  the  Pacific 
Basin  to  the  shallow,  runoff-freshened  reaches  of  the 
Yellow  Sea  and  East  China  Sea.  The  Kuroshio- 
Tsushima  Current  System  complexly  branches  and 
distributes  warm,  saline  waters  to  the  North  Pacific, 
Sea  of  Japan,  and  Yellow  Sea  after  transport  through 
and  mixing  within  the  East  China  Sea.  The  latter  two 
basins  have  vast  areas  of  shallow  water  of  <  lOOm 
depth,  strong  tidal  activity,  severe  atmospheric  forc¬ 
ing  events,  and  among  the  largest  freshwater  and 
fine-grained-sediment  fluxes  on  the  globe.  In  these 
shallow-water  regions  and  adjacent  shelf  edges,  hori¬ 
zontal  mixing  is  extremely  important  as  runoff  from 
coastal  regions  encounters  and  is  mixed  with  the 
equatorial  waters  of  the  Kuroshio-Tsushima  system. 
The  vigorous  winter-season  cooling  processes 
destratify  the  water  column,  forming  oceanic  fronts 
between  shallow  and  deep  waters  and  waters  of 
northern  versus  southern  origin.  These  effects 
facilitate  detection  and  tracking  of  various  water 
types  by  surface  temperature  using  satellite  infrared 
scanners  and  by  turbidity  (color)  using  satellite 
visual-range  sensors. 

Geographic  orientation  is  provided  by  Figure  3, 
illustrating  ocean  basins,  land  masses,  and  national 
boundaries.  Figure  4  is  a  daytime  infrared  image  of 
the  region  from  the  older  NOAA-5  satellite  of  8  April 
1978.  It  shows  the  patterns  of  earth-emitted  radia¬ 
tion,  with  highest  intensities  (dark)  over  the  Kuroshio 
and  sun-warmed  land  masses,  lowest  intensities 
(white)  over  the  cold  cloud  masses,  and  various  gray 
tones  showing  sea-surface  temperatures  with  well- 
developed  mesoscale  convolutions.  The  data  from 
this  older  satellite  system  had  a  spatial  resolution  of 
0.9  km  (at  nadir),  a  0.5°C  quantization  interval,  and 
an  analog  data  link  with  a  resulting  noise  level  of  1-2° 
RMS.  The  current  generation  of  NOAA  satellites 
(TIROS-N  series)  has  a  spatial  resolution  of  1.1  km, 
a  relatively  noise-free  digital  link,  and  a  0.125°C 


quantization  interval.  The  data  are  processed  here  to 
a  0.25  °C  resolution.  Both  systems  have  sun- 
synchronous  orbits  with  a  12-hour  repeat  period. 

Data  we  have  processed  from  one  of  the  new 
generation  of  satellites  (NOAA-6)  over  the  Korea 
Strait  are  shown  in  Figures  5-9.  These  are  parts  of 
two  overpasses  from  8  and  17  November  1980  ac¬ 
quired  in  computer-compatible  tape  format  from  Na¬ 
tional  Climatic  Center  and  analyzed  at  the  Remote 
Sensing  and  Image  Processing  Laboratory  of  Loui¬ 
siana  State  University.  Data  from  the  arrays  of 
measured  temperatures  acquired  by  the  satellite  are 
presented  in  three  formats:  imagery,  temperature- 
distance  transects,  and  spot  samples.  Spot  samples 
are  the  mean  values  of  temperature  within  the  tri¬ 
angular  areas  labeled  on  the  imagery.  The  images 
were  individually  enhanced  to  maximize  the  gray- 
shade  contrast  between  the  oceanic  features.  Imagery 
in  Figures  5,  7,  and  8  were  specially  processed  with 
very  cold  land  and  cloud  areas  electronically  black¬ 
ened  to  improve  photographic  reproduction.  Aside 
from  these  artificially  blackened  areas,  the  darker 
tones  represent  warmer  temperatures.  The  white  lines 
in  Figures  5  and  6  replicate  ship  track  lines  from  the 
Korea  Strait-East  China  Sea-Yellow  Sea  Remote 
Sensing  Experiment  (KEYRSEX  80),  which  was  con¬ 
ducted  in  November  and  December  1980.  During  this 
oceanographic  experiment  soundings  of  the  ocean 
and  the  atmosphere  were  obtained  for  correlati'  n 
with  satellite  imagery.  The  numbered  transec*  in 
Figure  8  do  not  indicate  ship  tracks,  but  are  pre¬ 
sented  to  show  surface  temperatures  and  the  strength 
of  the  gradients  observed  in  the  image. 

Overviews  of  both  scenes  in  Figures  5  and  6 
show  the  warm  waters  of  the  Tsushima  Current  ex¬ 
tending  from  the  south  to  a  strong  front  between 
Cheju  and  Tsushima  Islands.  During  this  experiment 
and  throughout  the  year  this  front  is  convoluted  into 
a  series  of  mesoscale  loops  and  meanders.  Near 
Tsushima  Islands  the  warm  waters  split  into  two 
branches.  The  main  branch,  the  Tsushima  Current 
proper,  flows  eastward  along  the  coast  of  Honshu, 
Japan,  toward  the  Tsugaru  Strait.  The  smaller  north¬ 
ward  branch  is  the  East  Korean  Warm  Current,  fol¬ 
lowing  a  broadly  meandering  path  parallel  to  the 
peninsula.  Radiation  temperature  differences  across 
the  front  on  both  days  are  shown  in  the  transects  of 
line  E  for  each.  Sea-surface  measurements  and  STD 
stations  acquired  along  line  E  on  17  November 
verified  the  presence  and  strength  of  the  front.  Table 
I  summarizes  a  comparison  of  ship  and  satellite  mea¬ 
surements.  17  November  was  a  relatively  warm  and 
humid  day  for  this  region  at  this  time  of  year,  with 
surface  air  temperatures  of  around  16  °C  and  total 
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Fiflure  4.  SOAA-5  saiellile  thermal  infrared  image  (10.5-12.5  mm.l  covering  the  Sea  of  Japan,  Yellow 
.Sea.  East  China  Sea,  and  a  small  pan  of  the  Sorth  Paci  fic.  8  April  1978.  A  surge  In  the  winter 
nnmsoon  has  brought  doad-free  polar  continental  air  out  over  the  entire  region.  Clouds 
(white)  occur  in  a  broken  mass  over  southern  Japan  and  in  a  broad  hand  '  the  southeast. 


precipitable  water  levels  of  1.5  g  cm  As  noted  in 
Table  I,  radiation  temperatures  of  both  water  ma.ssc.s 
are  lowered  by  several  degrees  from  actual  surface 


values,  but  the  temperature  of  the  warmest  waters 
was  affected  the  most.  The  lowest  temperatures  and 
salinities  were  found  by  surface  measurements  to  oc- 
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Figure  5.  \OAA-6  satellite  thermal  infrared  image  (channel  4.  /0.5-1  l.^mm/  of  the  south  coast  of 
Korea  and  Korea  Strait,  2325/,  S  \o\  emlyer  19S0.  Mean  temperature  of  the  contained  picture 
elements  is  labeled  next  to  the  triangular  sample  areas.  I  ettered  lines  are  locations  of  ship 
tracks  of  cruise  that  commenced  on  14  Soxemher.  Line  E  of  satellite  temperatures  versus 
distance  is  plotted  below.  Cold  land  and  cloud  tops  have  f/een  electronically  blackened  to  op¬ 
timize  the  display. 


cur  just  north  of  the  front  rather  than  in  the  coastal 
regions.  This  narrow  belt  of  cool,  low-salinity  water 
extends  from  the  East  China  Sea- Yellow  Sea  region 
west  of  Cheju  eastward  through  the  strait,  't  is  shelf 
water  known  to  discharge  eastward  between  the 
north  wall  of  the  Tsushima  Current  and  the  local 
coastal  water  from  late  summer  to  early  winter. 

Electronically  expanded  sections  of  these  images 
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are  presented  in  Figures  7  and  8.  The  narrows  of  the 
strait  by  the  Tsushima  Islands  from  8  November  are 
enlarged  to  show  mesoscale  details  with  greater  clar¬ 
ity.  The  radiation  temperature  patterns  and  values 
from  the  southern  Sea  of  Japan,  just  east  of  Korea, 
from  17  November  are  shown  in  Figures  8  and  9.  The 
extreme  enhancement  of  the  image  provides  water 
mass  outlines  and  temperature  stream  line  patterns. 
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The  land  masses  are  extremely  cold  at  this  hour  (7:30 
a.m.  local  sun  time)  and  so  produce  an  abrupt  tem¬ 
perature  drop  at  the  land-water  boundary.  Low 
cumulus  clouds  to  the  right  of  the  imagery  near 
Japan  generate  the  distinct  sawtoothed  pattern  seen 
in  transect  1  of  Figure  9.  From  this  perspective, 
Tsushima  Island  mass  acts  as  a  fixed  water  vane, 
deflecting  parcels  of  Tsushima  Warm  Current  waters 
from  the  northward  meanders  of  the  frontal  bound¬ 
ary  along  the  path  of  the  East  Korean  Warm  Cur¬ 


rent.  As  a  result,  alternating  warm  and  cold  water 
parcels  flow  north  along  the  east  coast  of  Korea. 

Conclusions 

A  practical  consequence  of  these  studies  is  the 
knowledge  that  usable  information  on  the  horizontal 
temperature  structure  of  the  upper  ocean  may  be  rea¬ 
sonably  expected  behind  outbreaks  of  polar  conti¬ 
nental  air.  Under  other  conditions,  however,  such  as 


6,  ,‘in  electronic  expansion  of  the  vicinity  of  Korea  Strait  front  the  previous  S  \ovemhcr  overpaw  Jhi\ 
suhuenc  includes  Korean  coast  to  teh  northwest,  Tsushima  islands  in  the  center,  and  cloud'^  pfio  oihcf 
Japanese  v^lands  to  the  southeast.  Sea-surface  temperature  pattern'^  \how  Tsushima  \\  ann  i  urrcni 
waters  interleaved  with  the  cold  waters  of  the  coastal  i  efiion. 
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Figure  7.  NOAA-6  sateUife  thermal  infrared  image  (channel  4,  10.5-1 1.5  mmi  of  the  Korea  Strait  and 
parts  of  adjacent  seas,  2326Z,  H  November  1980.  The  temperature-distance  plot  corresponds 
to  a  line  of  oceanographic  measurements  made  by  a  ship  on  the  same  date. 


in  humid  tropic  regions  or  deep  ocean  areas  where 
maritime  air  masses  predominate,  the  single-channel 
data  from  a  satellite  infrared  radiometer  (8-13  /»m 
band  of  the  spectrum)  has  little  chance  to  faithfully 
record  surface  temperature  patterns.  On  longer  time 
and  space  scales,  the  atmospheric  marine  boundary 
layer  (MBL)  is  strongly  influenced  by  the  sea  surface, 
its  major  source  of  heat  and  water  vapor.  Satellite 
observations  of  radiation  temperature  patterns  may 


follow  oceanic  fronts  in  humid  regions  by  virtue  of 
the  strong  influence  that  surface  conditions  have  on 
the  MBL.  As  sea-surface  features  move  or  change,  its 
thermal  counterpart  in  the  MBL  may  lag  behind,  or 
vice  versa.  We  would  judge  that  cloud-free  temper¬ 
ature  patterns  in  warm,  humid  oceanic  regions  would 
not  serve  as  an  operationally  reliable  link  to  the  up¬ 
per  ocean  because  of  the  strong  surface  atmospheric 
absorption  and  emission  effects  noted  above. 
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Figure  8.  NO  A  A -6  satellite  infrared  image  of  the  Korea  Strait  and  southern  Sea  of  Japan  eleclroniialh 
expanded  from  the  1 7  November  1980  overpass.  Flow  paths  of  the  Tsushima  and  East  Korean 
Harm  Currents  are  shown  in  the  sea-surface  temperature  patterns  of  the  southern  Sea  of 
Japan.  Temperature  transects  I  and  2  are  presented  in  Figure  9. 


However,  analogous  to  cold-air  outbreaks  are 
localized  mesoscale  winds  such  as  the  mistral  and 
land  breezes.  Although  they  are  smaller  in  spatial  and 
temporal  scales,  they  too  blow  from  land  to  sea.  The 
drier  and  cooler  air  characteristic  of  such  offshore 
winds  may  also  provide  “windows”  for  remote  sens¬ 
ing  of  coastal  and  continental  shelf  waters.  Certainly 
more  meteorological/oceanographic/remote  sensing 
team  studies  are  needed  to  explore  and  make  prac¬ 
tical  the  use  of  atmospheric  “windows”  on  the  sea. 

Avoidance  of  the  atmospheric  limiting  effects  is 
perhaps  possible  using  other  parts  of  the  spectrum  or 
using  a  multispectral  approach,  with  calculations 
from  spectral  attenuation  i  lerences  to  correct  each 
image  element.  The  mid-intrared  data  from  channel 
3  (3. 5-3. 9  ^im)  of  the  NOAA  satellites  is  much  less 
susceptible  to  atmospheric  effects,  and  preliminary 
studies  indicate  promise  for  this  sensor. ‘  '  It  is  usable 
at  night  only,  because  of  the  strong  solar  component 
in  this  part  of  the  spectrum.  At  this  time  it  appears 
that  true  all-weathi..  day-night  sea-surface  temper¬ 


ature  data  must  await  development  of  high-resolu¬ 
tion  multispectral  passive  microwave  radiometers.  In 
the  meantime  there  is  much  oceanography  to  be 
learned,  models  to  be  tested,  and  operational  support 
possible  under  the  outbreaks  of  polar  continental  air 
around  the  globe. 
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Figure  9.  NOAA-6  temperature-distance  Iransecis  I  and 2  from  Figure  8,  I '  November  1980.  The  north 
wall  of  the  Tsushima  Current  near  Japan  (upper)  and  the  wes'  wall  of  the  East  Korean  H  arm 
Current  (iower)  are  2.5  °C  radiation  temperature  steps.  Intersection  with  land  and  cloud  tops 
produces  the  lowest  temperature  values. 
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Field  Studj^  of 

“Anomalous  Cray-Shades*’ 


Satellites  have  been  used  for  many  years  to 
observe  clouds  and  cloud  systems  to  better 
understand  atmospheric  processes.  Clouds  are 
rather  easy  to  “see"  with  a  satellite  due  to  their 
high  reflectance  or  albedo.  Recently  introduced 
satellite  sensors  are  capable  of  detecting  variations 
of  very  low  albedo  in  the  atmosphere.  These 
variations  of  brightness  away  from  bright  cloudy 
regions  are  called  “anomalous  gray-shades"  and 
are  frequently  observed  in  visible  and  near- 
infrared  images  from  other  satellite  systems. 

A  team  of  investigators  under  the  direction  of 
Professor  Thomas  Yonder  Haar  at  Colorado  State 
University,  funded  by  the  Office  of  Naval 
Research’s  Selected  Research  Opportunity 
Program  has  been  studying  these  anomalous  gray- 
shades. 

To  investigate  these  features,  a  coordinated 
research  effort  was  conducted  between  19  Sep¬ 
tember  and  2  October  1980,  off  shore  of  Southern 
California.  Airborne  state-variable,  air  motion, 
particulate  and  radiometric  measurements  were 
obtained  to  determine  the  structure  of  gray-shade 
regions  identified  in  satellite  images.  Preliminary 
results  indicate  that  certain  of  these  features  are 
associated  with  the  marine  boundary  layer.  They 
seem  to  be  related  to  variations  in  size  and  con¬ 
centrations  of  large  particles  primarily  within  the 
boundary  layer.  (Figure  I  shows  a  set  of  features 
observed  with  the  NOAA-6  Advanced  Very  High 
Resolution  Radiometer  (AVHRR).  They  were 
observed  on  the  28th  of  September  in  the  wake  of 


Figurr  I.  Anomalous  Cray-Shade  Features  in  the  Wake  of 
Santa  Rosa  and  Santa  Cruz  Islands  off  of  Santa 
Barbara,  California  as  Observed  by  the  NOAA-6 
Advanced  Very  High  Resolution  Radiometer 
(A  VHRR)  on  28  September  1980. 

Santa  Rosa  and  Santa  Cruz  Islands  just  off  of 
Santa  Barbara,  California.  The  features  are 
aligned  with  the  boundary  layer  wind  and 
although  some  of  the  brightness  variations  can  be 
attributed  to  surface  roughness  variations,  most 
of  the  signal  is  due  to  atmospheric  reflection.) 
Anomalous  gray-shade  regions,  when  fully  un¬ 
derstood,  may  provide  great  observational  insight 
into  stratocumulus  formation  and  dissipation, 
humidity  variations,  visibility  determination  and 
forecasting  as  well  as  boundary  layer  dynamics 
and  air-sea  interactions.  ■ 
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Surface  Wave  modulation 
Patterns  in  Radar  Images 


by 


Erik  nflolb-Christer^n 
fOossachusetts  Institute  of  Technology 


Images  formed  from  radar  sea  surface  scatter 
return  to  aircraft  and  satellite  platforms  show 
coherent  patterns  that  can  be  associated  with  current 
shear,  eddies,  internal  waves,  and  patterns  on  the  sea 
surface  and  bottom  topography.  Seasat  SAR  (Syn¬ 
thetic  /Iperture  /?adar)  data  show  a  wealth  of  de¬ 
tailed  patterns.  An  example  of  this  in  Figure  1,  shows 
a  digitally-processed  image  of  Nantucket  Island, 
Massachusetts  and  its  surroundings.  The  magnified 
part  in  Figure  2  indicates  a  velocity  shear  at  the  edge 
of  a  current  of  less  than  one  knot  running  to  the 
WNW,  south  of  the  Shoals.  One  cannot  explain  these 
changes  in  the  wave  packet  structure  of  sea  surface 
waves  in  terms  of  the  results  of  the  customary  ap¬ 
proximation  in  terms  of  waves  of  infinitesimal  ampli¬ 
tude  or  linear  theory.  The  patterns  in  wave  amplitude 
and  wavelength  are  due  to  effects  of  finite  wave  am¬ 
plitude,  where  nonlinear  phenomena  become  impor¬ 
tant.  The  tendency  for  waves  to  undulate  in  ampli¬ 
tude  and  frequency  is  typical  of  finite  amplitude 
waves.  We  refer  to  these  as  amplitude  (AM)  and  fre¬ 
quency  (FM)  modulations,  and,  since  they  are  due  to 
nonlinear  dynamics,  we  call  them  nonlinear  wave 
modulations. 

I  shall  attempt  to  describe  in  this  article  some 
features  of  modulational  dynamics  for  surface 
waves,  and  how  the  detection  process  used  for  SAR 
and  Side  Looking  Airborne  Radar  data  rely  on  inter¬ 
modulations  of  the  scattering  wavelets  and  the  under¬ 


lying  dominant  wave  field.  Then  1  will  go  on  to  sug¬ 
gest  that  the  same  kind  of  forced  modulation  and  in¬ 
termodulation  changes  the  patterns  in  the  dominant 
wavefield  to  make  eddies,  shear  layers  and  other 
coherent  motions  visible.  First  consider  modulations 
of  the  surface  wavefield. 

Surface  Wave  Modulations 

Both  surfers  and  boatmen  know  that  swell  and 
wind  waves  vary  in  amplitude  and  arrive  in  packets 
or  groups.  A  wavefield  tends  to  develop  variations  in 
wave  amplitude,  a,  because,  at  a  given  wave  steep¬ 
ness,  a  dip  in  amplitude  (a,,^  >  0)  tends  to  move 
faster  than  a  positive  bump  in  wave  amplitude  (a,^,^  < 
0).  So,  dips  and  bumps  will  tend  to  separate  rather 
than  merge  and  cancel.  This  instability  process  makes 
disturbances  in  wave  amplitude  grow  initially,  when 
a  uniform  amplitude  wavefield  is  disturbed.  But  near 
the  bumps,  there  tends  to  be  larger  wave  steepness 
ak  =  27ra/X,  defined  in  terms  of  amplitude,  a,  and 
wavelength,  \.  Steeper  waves  tend  to  propagate 
faster,  so  there  will  be  a  tendency  for  the  front  of 
wave  groups  to  steepen  further. 

These  effects  of  amplitude  variation  and  steep¬ 
ness  on  the  phase  speed  of  finite  amplitude  periodic 
waves  (Stokes  waves)  are  contained  in  the  dispersion 
relation,  relating  phase  speed  c,  wave  number  k  = 
2ir/X,  and  wave  steepness  ak  =  e 


s 


r 


gk(l+a'k'  + 

=  gk(l  +  +  a„/(8ek)l. 

Eventually,  in  the  interplay  between  steepness  and 
bumpiness  effects,  an  equilibrium  situation  is  estab¬ 
lished,  so  the  mature  wave  field  consists  of  constantly 
changing  groups  of  waves,  maintaining  an  average 
steepness  and  wave  amplitude.  This  is  the  typical  sea 
surface  wave  field  under  weakly  generating  wind 
wave  conditions,  with  nearl>  constant  average  prop¬ 
erties,  but  always  changing  locally  in  terms  of  group 
lengths  and  amplitude.  Now  and  then  accidents 
happen,  in  the  sense  that  a  further  instability  occurs 
where  subharmonic  transitions  take  place,  and  a  crest 
is  lost.  Due  to  such  accidents,  the  average  wavelength 
of  a  wavefield  will  tend  to  increase  and  the  wave 
steepness  decrease  until  there  is  an  open  ocean  swell 
wavefield  of  modest  steepness.  The  description  I 
have  given  is  based  on  analytical  results  for  nonlinear 
waves,  but  not  all  the  detailed  steps  have  been  veri¬ 
fied  as  occurring  in  the  order  described,  so,  in  that 
sense,  my  description  is  somewhat  speculative. 


Figure  /.  Digitally  processed  image  of  Nantucket  Island  and  the 
surrounding  waters  obtained  from  SEA  SAT  SAR  data. 


Linear  and  Nonlinear  Wave  Dynamics 

Linear  theory  is  excellent  for  describing  wave 
development  over  limited  distances  and  times;  for 
surface  waves,  linear  approximations  are  good  for 
propagation  distances,  L,  of  much  less  than  X/e',  for 
larger  distances  nonlinear  effects  become  more  im¬ 
portant.  For  100  m  wavelength  swell  of  amplitude  1 .5 
m  or  nearly  5  ft.  crest  to  trough  height,  e  =  0.1  and 
the  propagation  distance  is  10  km.  For  small  wavelets 
of  10  cm  wavelength  and  0.5  cm  amplitude  e  =  0.3, 
the  propagation  distance  is  only  one  meter. 

For  propagation  over  longer  distances,  the  cu¬ 
mulative  effects  of  nonlinear  terms  in  the  equations 
of  motion  become  important,  and  one  has  to  develop 
the  dynamics  to  third  order  in  wave  steepness  and  ob¬ 
tain  the  Davey-Stewartson  equation’  in  wave- 
complex  amplitude  A(x,y,t), 


'(A,  +  CgA,) 

“  ^2  (■Nx  ~  (2) 

-  -2  |A|>  A  =  0. 


Figure  2.  Magnified  portion  of  Figure  /.,  showing  the  shear 
instability  of  a  current  edge  shear  layer. 
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The  water  surface  deflection  is 

z  =  ffx.y.t)  =  Re|A(x,y,t)  e**''*  -  |.  (3) 

Here,  k  and  w  are  the  average  wavenumber  k  =  2t/X 
and  average  frequency  w  =  2ir/T,  T  is  the  wave 
period. 

The  Davey-Stewartson  equation  gives  a  valid  ap¬ 
proximation  over  distances  of  propagation  much 
smaller  than  Li  =  X/e*.  (Lj  is  not  the  same  as  the  L 
for  linear  approximation.  It  gives  a  “good”  descrip¬ 
tion  for  propagation  time  scales  of  order  T/e^  and  for 
group  lengths  of  order  X/r.  For  the  swell  example, 
X  =  100  m,  t  =  0. 1 ,  Lj  =  100  km.  For  the  wavelets 
of  X  =  10  cm  and  «  =  0.3,  Ls  =  30  m. 

The  Davey-Stewartson  equation  (2)  has  exact 
solutions^  in  the  form  of  periodic  and  solitary  distur¬ 
bances.  In  two  dimensions,  this  equation  reduces  to 
the  cubic  Schrodinger  equation*  *,  first  used  in  wave 
mechanics  to  solve  for  the  electron  as  a  wave  packet. 
Also  when  applied  to  sea-surface  wave  modulations, 
the  solutions  to  the  cubic  Schrodinger  equation  have 
quantum  properties,  limiting  the  energy  content  of  a 
single  wave  packet  for  a  given  wavelength.  The  solu¬ 
tions  also  possess  certain  conservation  properties,  so 
that  wavepackets  that  collide  will  survive  collisions. 
This  suggests  that  wave  packets  can  have  long  half 
life,  and  that  wavepacket  descriptions  of  a  sea- 
surface  wave  field  may  be  useful  for  some  purposes.* 
But,  as  will  be  argued  further  along,  the  modulation 
field  is  also  very  sensitive  to  disturbances,  which  can 
alter  the  wave-packet  structure  of  a  wave  field  en¬ 
countering  current,  wind  or  topography. 

Forced  Modulation 

The  modulation  field,  as  described  by  equation 
(2),  will,  when  otherwise  undisturbed,  reach  an 
equilibrium  state  where  the  linear  and  nonlinear 
terms  balance.  A  current  shear  dv/dx  =  vk,  where 
the  mean  velocity  changes  by  an  amount  v  in  a  wave¬ 
length  X,  will  introduce  additional  terms  in  the 
Davey-Stewartson  equation.  One  term  comes  from 
the  energy-flux  divergence  term  in  the  energy  equa¬ 
tion*,  and  introduces  a  term  (vA),  »  (vAk)  in  equa¬ 
tion  (2).  There  are  additional  terms  due  to  radiation 
stresses',  but  they  are  of  the  same  order-of- 
magnitude. 

Velocity  shear  will  be  as  important  as  the  non¬ 
linear  term  in  equation  (2)  when 

Avk  *  — -  |A|*  A, 

2 


setting 

Cg  =  and  k'  (Aj*  =  (\ 
2k 

one  finds  that 


is  the  criterion  for  current  shear  to  affect  modula¬ 
tions.  This  shows  that  small  current  shear  can  be  im¬ 
portant. 

For  an  example  of  a  100  m  wavelength  swell  of 
steepness  e  ~  0. 1 ,  the  group  velocity  is  Cg  =  6  m/s 
and  the  current  variation  required  to  disturb  the 
modulation  field  is  v  =  t’Cg  s  6  cm/sec  =  .12  knots. 
For  the  6  cm  waves  of  steepness  e  =  0.3,  the  current- 
change  needed  is  of  order  v  =  1.5  cm/sec.  For  less 
steepness,  with  r  =  0.1  and  X  =  6  cm,  we  find  the 
current-variation  needed  to  upset  the  modulation 
field  to  be  of  order  v  s  1 .5  mm/s. 

This  demonstrates  how  very  small  coherent  dis¬ 
turbances  in  the  current  field  can  cause  disturbances 
in  wave  modulation  patterns  in  wavefields  of  moder¬ 
ate  steepness,  and  if  the  disturbances  are  coherent, 
the  modulational  disturbances  will  also  be  coherent. 
TTie  steeper  the  waves,  the  larger  the  disturbance 
needed  to  affect  the  modulation  field.  But  ship  wakes 
and  weak  eddies  will  show  up  in  the  modulation  field 
of  short  wavelets  of  X  =  6  cm  under  a  wide  range  of 
velocity  disturbances  down  to  a  few  cm/s. 

Disturbances  in  wave  propagation  speed  due  to 
bottom  topography  will  affect  the  modulation  field 
the  same  way;  for  shoaling  from  very  deep  to  a  depth 
of  20  m  will  change  the  group  velocity  of  100  m 
waves  from  6,25  m/sec  to  5.32,  so  that  the  change  is 
V  =  0.93  m/sec,  sufficient  to  change  the  modulation 
field.  This  suggests  the  reason  why  one  can  see  the  lo¬ 
cation  of  the  continental  shelf  edge  in  some  SAR  im¬ 
ages,  and  why  topography  shows  up  so  very  clearly. 

So  far,  we  have  discussed  how  waves,  as  they 
propagate,  eventually  accumulate  effects  of  nonlin¬ 
earity  until  they  add  up  enough  to  change  the  wave- 
field  by  introducing  modulation  in  phase  and  ampli¬ 
tude.  The  steeper  the  waves  are,  the  quicker  they 
show  the  effects  of  nonlinear  processes.  But  through 
all  this,  the  linear  approximation  is  a  good  one  for 
small  time  and  distance  of  propagation.  So  observa¬ 
tions  of  wavefields  over  short  distances  and  times  will 
tend  to  be  dominated  by  linear  effects  and  reveal  lit¬ 
tle  of  the  nonlinear  dynamics.  But  for  the  much 
larger  scales  of  observation  made  possible  by  remote 
sensing,  the  nonlinearities  have  room  to  show  up. 
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and  give  disturbances  their  chance  to  modify  wave 
modulation  fields. 

Intermodulation 

The  current  that  disturbs  the  wave  modulation 
field  may  be  a  fluctuating  current  due  to  other  waves. 
A  field  of  \  =  10  cm  waves  riding  on  waves  of  100  m 
wavelength  of  3  m  amplitude  (t  =  0.10)  will  be  sub¬ 
ject  to  current  fluctuations  of  ±  1.2  m/sec.  Accoiu- 
ing  to  our  calculations,  this  current  is  sufficient  to 
dominate  the  modulation  field  of  the  10  cm  waves. 
The  orbital  velocity  field  is  convergent-divergent,  so 
that  it  plays  accordion  with  the  wavelet  field,  induc¬ 
ing  large  steepness  and  leading  to  fourth-order  ef¬ 
fects  becoming  important,  the  wavelet  field  may 
undergo  subharmonic  transitions  and  even  wave¬ 
breaking  at  certain  phases  of  the  underlying  domi¬ 
nant  sea. 

The  mechanism  of  forced  modulation  of  the 
wavelets  by  the  dominant  sea  is  the  basis  for  the  SAR 
detection  of  sea  surface  waves*  ’  The  SAR  signal 
is  scattered  from  wavelets  and  the  intensity  and  Dop¬ 
pler  shift  of  the  reflected  wave  is  interpreted  as  a 
modulation  signal  due  to  the  dominant  wave  and, 
through  a  suitable  algorithm,  translated  as  sea  sur¬ 
face  shape  due  to  the  underlying  waves.  This  inter¬ 
pretation  algorithm  for  SAR  and  other  radar  scatter¬ 
ing  imagery  is,  of  course,  a  crucial  element  in  data 
interpretation.  Because  of  the  complexity  of  the  in¬ 
termodulation  processes,  one  has  to  develop  better 
algorithms  by  simultaneous  observation  of  both 
wavelets  and  underlying  waves  over  a  range  of  wind 
and  wave  conditions.  This  should  preferably  be  done 
without  the  need  for  flying  radars  on  aircraft  or 
satellites,  since  the  cost  of  sufficiently  extensive 
overflights  would  be  prohibitive.  I  have  therefore 
designed  a  scattering-wave  gauge  which  simulates 
radar  backscatter. 

The  “Harp”  Wave  Gauge 

A  student,  J.  Menoher,  and  I  have  built  and 
tested  a  simple  wave  gauge  whose  output  simulates 
radar  scatter  return.  The  wave  gauge  consists  of  a 
single  capacittince  wire,  strung  on  a  frame  with 
parallel  vertical  segments  a  distance  f  apart,  as  shown 
in  Figure  3.  The  capacitance  of  the  wire  varies  pro¬ 
portionally  to  the  sum  of  the  water  level  displace¬ 
ments  at  all  the  wires.  A  capacitance  detector  circuit 
will  then  provide  a  wave  gauge  output  proportional 
to  the  sum  of  displacements.  For  very  long  wave¬ 
lengths,  compared  to  the  length  of  the  array,  the 
gauge  will  work  like  a  single  point  wave  gauge.  For 


Figures.  **Harp”  wave  gauge,  with  a  wire  strung  in  parallel  runs. 
Resonant  wave  contour  shown. 


waves  of  wavelength  equal  to  the  wire  spacing,  the 
capacitance  variation  of  all  the  wires  will  be  in  phase, 
and  the  output  will  be  a  maximum.  For  wavelengths 
twice  the  wire  spacing,  the  output  will  be  zero.  The 
output  amplitude  will  vary  with  wavelength  as  shown 
in  Figure  4.  The  wave  gauge  functions  as  a  resonant 
scattering  detector.  It  can  also  be  built  for  less  than 
the  cost  of  a  satellite;  its  capital  cost  will  be  com¬ 
parable  to  an  aircraft  flight  hour. 

A  first  test  of  a  ten-wire  wave  “harp”  gave  the 
results  shown  in  Figure  5.  The  upper  trace  is  the 
“harp”  output  and  the  lower  one  is  the  output  of  a 
single-wire  capacitance  wave  gauge  placed  next  to  it 
in  a  wave  tank.  The  frequency  of  the  wave-generating 
paddle  increases  slowly,  and  the  resonance  peak  is 
evident.  But  it  is  also  interesting  to  note  that  there  are 
bursts  of  high-frequency  output  for  incoming  waves 
of  frequencies  of  1/2  and  1/3  of  the  frequency  of  the 
resonant  wave. 

This  means  that  a  scattering  detector  also  detects 
higher  harmonics  of  Stokes  waves.  The  same  must  be 
the  case  for  a  Bragg  scatter-return  of  radar  signals. 
This  result  illustrates  the  complexity  of  interpretation 
of  scatter  returns,  and  indicates  the  need  for  careful 
observations  to  establish  relationships  between  scat¬ 
tering  waves,  wind  and  dominant  sea.  Data  process¬ 
ing  of  data  from  such  an  experiment  or  from  SAR 
data,  for  that  matter,  requires  the  use  of  processing 
methods  that  produce  information  about  modula¬ 
tions  in  phase  and  amplitude  and  intermodulation;. 

Eduction  of  Modulations  from  Signals 

What  we  need  is  information  about  the  complex 
modulation  variable  A  in  equation  (3).  One  method, 
taken  from  communications  engineering,  uses 
Hilbert  Transforms"  and  has  been  used  for  surface 
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Figure  5.  Output  of  "harp”  wave  gauge  (upper  trace)  compared  to  signal  from  a  single  wave  gauge 
dower  trace)  as  wave  frequency  is  slowly  increased  showing  higher  harmonic  resonance 
lAAB)  and  main  resonance  peak  (C).  (Result  obtained  by  J.  Menoher). 
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wave  processing  by  my  colleague  W.  K.  Melville.  One 
of  our  students,  G.  Sahar,''  has  used  the  method  to 
process  laboratory  wind  wave  data. 

Figure  6  shows  a  wave  gauge  output,  with  the 
amplitude  !A|  superimposed,  and,  in  the  lower 
trace,  the  phase  ©  of  A  =  Re*®;  the  linear  trend  due 
to  the  carrier  wave  txp[i(kw-a)t)]  has  been  taken  out. 
Once  one  obtains  A,  one  can  relate  it  to  the  other 
scale  processes  that  drive  or  interact  with  the  modula¬ 
tions.  The  interpretation  of  SAR  return  by  scattering 
wavelets  and  relating  it  to  the  dominant  sea  is  a  dif¬ 
ficult  problem  and  is  similar  to  the  challenges  of 
relating  wave  modulation  to  internal  waves  and  other 
large  scale  processes.  They  present  similar  data  proc¬ 
essing  problems. 

In  learning  to  interpret  data  in  terms  of  such  in- 
termodulational  relationships,  we  need  to  do  some 
controlled  experiments,  guided  by  existing  work  on 
nonlinear  waves  and  intermodulation  and  with  an 
awareness  of  the  limitations  of  simple  and  convenient 
approximations. 

Applications  to  Oceanography 

The  real  advantage  of  active  electromagnetic 
sensors  in  oceanography  is  their  insensitivity  to 
clouds  and  fog.  One  does  not  realize  the  prevalence 
of  cloud  and  fog  patches  before  one  starts  looking 
for  visible  and  infrared  satellite  data  for  the  sea  sur¬ 
face.  For  some  locations,  such  as  the  Drake  Passage, 
for  example,  one  may  find  one  day  a  year  with  a 
good  view.  Other  areas  are  nearly  as  frustrating; 
climatological  records  of  sunshine,  cloud-cover  and 
fog  are  readily  available  to  demonstrate  this.  Also, 
interesting  events,  such  as  hurricanes,  gales  and  in¬ 
tense  air-sea  interaction  do  not  occur  on  sunny  days. 
So,  without  sensors  that  see  through  clouds,  we  can¬ 
not  obtain  the  synoptic  view  of  patterns  that  are 
uniquely  available  from  remote  sensing  devices. 

Good  distance-resolution,  another  crucial  fea¬ 
ture  of  SAR  and  SLAR  data,  derives  from  the  fact 
that  distance  is  measured  in  terms  of  time  for  a  signal 
to  return  rather  than  through  optical  resolution  of 
geometry.  A  spatial  resolution  better  than  fifty 
meters  makes  it  possible  to  see  the  dominant  surface 
waves  and  swell,  internal  waves  and  current  shear. 
Because  the  dominant  waves  can  be  resolved,  it  is 
also  possible  to  educe  from  the  data  the  amplitude 
and  modulation  field  of  those  waves,  and  to  interpret 
the  modulation  patterns  in  terms  of  larger  scale 
phenomena.  It  is  worth  emphasizing  that,  because 
the  modulation  field  is  governed  by  nonlinear  proc¬ 
esses,  one  cannot  use  linear  filtering  theory  to 
calculate  resolution.  The  resolution  of  an  image 


needed  to  discover  a  100  km  eddy  from  the  surface 
wave  modulation  field  must  be  sufficient  to  resolve 
individual  waves. 

As  an  alternative,  one  can  build  a  modulation 
eduction  operator  into  the  processing  scheme.  One 
reason  for  the  clear  view  of  amplitude  modulations 
afforded  by  existing  SAR  imagery  is  that  somehow 
such  an  eduction  operator  is  present  in  the  data  proc¬ 
essing  chain,  possibly  by  inadvertence  or  because  of 
the  skewness  of  modulations  of  scatterers  on  surface 
waves . 

I  will  next  give  a  few  examples  of  data  use  in 
oceanography. 

Density  Structure  from  Internal 
Wave  Patterns 

When  and  if  SAR  and  similar  active  sensor  sea 
surface  data  become  available,  it  will  be  possible  to 
put  into  routine  use  a  method  for  estimating  vertical 
density  structure  that  a  student,  A.  da.  S.  Mascaren- 
has,  Jr.,  and  1  developed  a  few  years  ago.  We  used 
Landsat  images  that  showed  patterns  that  could  be 
identified  with  internal  waves  propagating  on  the 
seasonal  thermocline.  These  internal  wave  packets 
occur  in  the  New  England  continental  shelf  region. 
The  waves  are  generated  by  tidal  flow  over  topog¬ 
raphy,  and  the  packets  are  generated  once  per  Ms 
tidal  period,  which  is  approximately  12.4  hours.  One 
sees  packets  generated  several  tidal  periods  apart  in 
many  Landsat  images.  The  spacing  between  packets 
is  the  distance  travelled  in  12.4  hours.  This  gives  the 
propagation  speed,  c,  of  the  wave  packets. 

The  linear  approximation  for  internal  waves  on 
a  density  jump  Aq/q  at  a  depth  h  from  the  surface  in 
otherwise  deep  water  gives 


Since  we  have  found  c  from  the  satellite  image,  one 
can  calculate  h  .  If  one  knows  Aq/q  from  other 
sources,  one  can  find  h.  We  argued  that  the  density 
jump  on  the  seasonal  thermocline  was  due  to  the  tem¬ 
perature  difference,  so  that  Aq/q  =  -  aAT.  We 

then  said  that  hAT  was  a  measure  of  heat  storage  in 
the  upper  layer,  and  calculated  its  distribution  for  the 
months  of  July  and  August'*.  (Figure  7). 

We  were  constrained  by  limited  data  availabil¬ 
ity,  caused  by  fog  and  the  Landsat  data  acquisition 
schedule.  This  method  may  well  prove  applicable  to 
internal  waves  on  the  main  thermocline,  since  one 
should  expect  such  waves  to  be  generated  all  along 
the  continental  edges  by  the  tides.  We  may  need  bet- 
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Figure  6.  Wave  gauge  output  from  laboratory  wind  wave  ex¬ 
periment,  with  amplitude  modulation  superimposed 
and  phase  modulation  shown  in  lower  trace.  (Results 
obtained  by  G.  Sahar). 


Figure  7. 

Heat  storage  in  upper  ocean  derived /tom  tntemal  wave 
packet  pacings.  MoRoChrtetensen  (I9tlh- 


ter  methods  of  processing  radar  data  to  see  the  sur¬ 
face  modulation  patterns  generated  by  these  internal 
waves.  The  distance  between  wave  packets  will  give 
the  speed,  and  one  can  again  observe  h  .  Another 
phenomena  that  can  be  excited  by  the  tides  is  shear- 
layer  instability. 

Estimating  Current  Speed  from 
Satellite  Data 

Another  example  of  interpretation  of  patterns 
visible  in  radar  images  is  given  by  regarding  the  un¬ 
dulations  of  a  current  edge  shear  layer  as  a  tidally 
forced  hydrodynamic  instability.*’ 

Where  a  current  passes  a  prominent  topography 
feature  and  speeds  up  locally,  oscillations  and  tur¬ 
bulence  tend  to  decrease  in  comparison  to  the  orga¬ 
nized  motion;  in  part,  this  is  due  to  the  increase  in 
mean  velocity,  without  a  corresponding  increase  in 
oscillations  and  turbulence.  Also,  since  the  current 
touches  a  boundary,  there  will  be  little  oscillation 
normal  to  the  boundary.  Further  downstream,  where 
the  current  loses  the  support  of  the  boundary  and  the 
current  edge  becomes  a  free  shear  layer,  the  flow 
becomes  unstable.  An  unstable  flow  tends  to  amplify 
fluctuations  due  to  external  disturbances,  so  that  the 
shear  layer  develops  lateral  fluctuations  in  the  form 
of  travelling  waves  of  increasing  amplitude.  This  is 
clearly  visible  in  TIROS-N  infrared  images  of  the 
Gulf  Stream  near  Cape  Hatteras,  North  Carolina, 
shown  in  Figure  8.  Another  current  edge  can  be  seen 
in  the  Seasat  SAR  image  of  the  Nantucket  Shoals, 
where  the  edge  of  the  WNW  current  south  of  the 
Shoals,  reported  in  the  U.S.  Coastal  Pilot,  forms  an 
undulating  shear  layer  edge,  shown  in  Figure  2. 

The  excitation  for  shear  flow  instability  will  be 
velocity  fluctuations  caused  by  incoming  internal 
waves,  storms  and  tides.  The  most  reliable  excitation 
is  the  tide,  and,  in  the  examples  shown,  the  M,  tide. 

We  can  calculate  the  speed  of  the  shear  layer 
waves  by  measuring  the  wave  length  in  the  satellite 
image  and  dividing  by  the  period  of  the  Mj  tide.  For 
the  Gulf  Stream  just  beyond  Cape  Hatteras  we  find 
c  «  1  m/sec.  Next,  we  use  the  fact  that  shear  layer 
waves  tend  to  move  at  a  speed  c  =  (U,  +  UO/2,  the 
average  of  the  current  velocity  on  the  two  sides.  We 
now  assume  that  the  current  velocity  on  one  side  of 
the  shear  layer  is  zero:  Uj  =  0;  one  then  finds 

II  ~  jr  -  2  (wavelength) 

^  (M-period) 


For  the  Gulf  Stream,  we  found  U,  =  2  m/sec,  very 
close  to  the  reported  value. 

In  order  to  verify  such  methods  of  current  esti¬ 
mation,  one  needs  a  series  of  sea-surface  images 
spaced  a  day  apart  or  so,  so  that  one  can  follow  a 
shear  layer  wave.  But  regions  where  a  cold  or  warm 
current  passes  a  headland  tend  to  be  active  mete¬ 
orologically,  breeders  of  storms,  regions  of  cumulus 
convection  or  just  foggy.  It  is  hard  to  find  infrared 
radiometer  images  of  Hatteras  or  the  Nantucket 
Shoals  free  of  clouds.  The  same  goes  for  most  other 
major  currents  like  the  Kuroshio,  the  Antarctic  Cir¬ 
cumpolar,  the  East  Greenland,  and  the  Somali  Cur¬ 
rent  as  well  as  the  many  minor  currents  we  could 
analyze  this  way. 

Again,  active  sensor  data,  such  as  SAR  or  SLAR 
data,  would  be  very  important  for  estimating  cur¬ 
rents,  analyzing  the  growth  of  fluctuations  and 
meanders  as  the  first  stage  in  formulating  predictive 
models  of  such  undulations. 

Ship  Wakes  as  Diagnostic  Tools 

Ship  wakes  show  up  in  high  resolution  images 
such  as  those  produced  by  SAR.  The  deformation  of 
a  ship  wake  by  current  shear  can  be  used  to  deduce 
the  shear,  provided  the  ship  keeps  steady  course  and 
speed,  and  one  knows  the  speed.  One  can  find  the 
speed  from  the  wavelength  of  the  stern  wave.  As 
anybody  who  has  spent  time  leaning  on  the  stern  rail 
of  a  ship  knows,  the  stern  wave  moves  at  the  speed  of 
the  ship;  this  means  that  the  stern  wave  propagates  at 
the  speed  of  the  ship,  v.  If  X  is  the  wavelength  of  the 
stern  wave,  one  has  the  relation 


V  s  s  (gX/2ir). 

Combining  the  analysis  of  shear  instability  with  use 
of  ship  wakes  as  a  diagnostic  can  aid  one  in  mapping 
out  current  fields  and  fluctuations  in  a  region.  This 
requires  high-resolution  sea-surface  data  without 
cloud  or  fog  clutter. 

The  Data  Volume  Problem 

High  resolution  data  represent  high  data  density 
per  unit  area  of  the  sea  surface.  SAR  data  with  a 
resolution  of  15  x  15  m’  contains  approximately  30 
megabytes  for  a  100  km“  data  field.  Transmission, 
storage  and  processing  of  such  large  data  volumes  are 
just  now  becoming  routinely  possible  at  manageable 
cost.  The  technologies  of  multiprocessors  and  mass 
storage  devices  are  advancing  so  rapidly  that  within  a 
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Jviia  haiKllnii!  pu>hU■lll^  and  costs  will  have 
sliaiiccd  lot  ilic  f'ciict,  and  wc  tinisi  trv  (o  aii(icit>a(c 
icchiiolocical  advance  as  best  we  can.  and  cness  at 
wheie  Inline  liolilenecks  will  i>ccur. 

Data  areliives  and  data  disiribiiiion  lo  research 
woikcis  ma>  end  up  on  hii;h-densii\  disk  media, 
while  operaiional  users  will  eel  ilieir  data  via  com- 
mumcaiioiis  saielliie. 


fhe  nonlinear  meihods  lor  lime  series  are  noi 
always  easy  lo  eenerali/e.  One  need  onl>  realize  ihai 
while  Ihe  order  in  which  one  applies  linear  processors 
ii>  a  signal  may  noi  alteci  ihe  resuli,  ihe  order  ot  ap- 
plieaiioii  of  nonlinear  processtrrs  can  make  a  large 
diflerence.  i  or  example 

sine  nil)  ^  \  sinuli)  and  log  nil)  t  (logu  i)e 


Improvemenl  in  Quantitative  Interpretation 

of  Data  Conclusion 


A  crucial  elemem  in  interpretaiion  ol'  sea- 
surlace  radar-scaiier  daia  are  nonlinear  wave  proc¬ 
esses.  We  need  to  know  such  processes  belter,  and  we 
need  to  he  able  to  analy/e  the  data  Held  to  take  apart 
the  results  of  nonlinear  interactions.  We  can  learn 
much  about  nonlinear  waves  from  laboratory  and 
field  experiments  and  much  of  the  knowledge  we 
need  can  be  obtained  without  havitig  to  Hy  satellites 
or  aircraft. 

The  analysis  methods  we  now  have  work  with 
time  series  and  joint  time  series,  that  is  for  funciions 
of  a  single  independent  variable.  We  need  to  extend 
such  methods  to  two-dimensional  and  three-dimen¬ 
sional  data  fields,  such  as  those  we  obtain  by  remtiie 
sensing. 


Radar  scatter  data,  are  now  useful  in  operations 
and  research.  Improvemenl  in  quantitative  inter 
pretation  of  the  data  requires  research  on  wave 
dynamics  and  data  processing,  lechnology  of  data 
transmission  and  storage  will  keep  improving  and 
will  be  available  to  handle  remote-sensing  data  when 
needed. 

Patterns  in  remote-sensing  images  can  be  inter¬ 
preted  quantitatively  to  produce  information  on  cur¬ 
rent  speed  and  density  structure.  Tor  such  use.  the 
all-weather  capability  and  the  good  spatial  resolution 
of  active  sensor  data  will  be  especially  valuable.  ■ 
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Introduction 

T'he  evolution  of  computer  technology  has  pro¬ 
vided  scientists  and  engineers  with  the  capability 
to  process  information  efficiently  in  the  form  of  im¬ 
ages  or  pictures.  This  has  resulted  in  a  great  deal  of 
research  interest  in  image  processing  as  exemplified 
by  the  May,  1981  special  issue  of  the  Proceedings  of 
the  IEEE'  devoted  to  papers  on  image  processing. 
The  collection  of  papers  emphasized  two  aspects  of 
image  processing:  (1)  the  processing  of  image  data  to 
achieve  desired  improvements,  and  (2)  image  under¬ 
standing  which  is  concerned  with  the  way  humans 
visualize  and  utilize  image  information. 

In  the  meteorological  satellite  community,  scien¬ 
tists  now  routinely  use  computers  to  process  image 
data  to  improve  image  quality  and  image  understand¬ 
ing.  For  example,  an  algorithm  to  remove  white 
background  due  to  snow  provides  a  better  quality 
picture  for  cloud  studies.  Another  example  is  the 
assignment  of  separate  colors  to  cloud  heights  which 
allows  the  user  to  differentiate  clouds  quickly  at  dif¬ 
ferent  levels  and  thus  provides  improved  image  under¬ 
standing.  The  two  examples  deal  with  flat  images  as 
viewed  from  the  satellite.  Currently,  the  majority  of 
satellite  imagery  is  viewed  in  this  context  and  the  user 
must  visualize  the  three-dimensional  structure  of  the 
atmosphere  in  his  mind. 

Given  a  flat  satellite  image,  contouring  algo¬ 
rithms  can  be  employed  to  operate  on  the  visual  or 
infrared  data  to  generate  three-dimensional  cloud  im¬ 
ages.  Another  alternative  is  to  use  polygons  to  ap¬ 
proximate  a  three-dimensional  cloud  surface.  These 
methods  result  in  three-dimensional  surfaces  that 
more  realistically  depict  what  a  pilot  will  see  from  an 
airplane  or  what  an  observer  will  see  from  the 
ground.  The  result  is  improved  image  utilization.  We 
will  define  image  utilization  as  a  transformation  that 
results  in  a  product  that  directly  depicts  a  desired 
phenomenon. 
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This  paper  describes  the  use  of  satellite  pictures 
and  computer-driven  imagery  to  represent  clouds  as 
three-dimensional  objects  in  a  viewer’s  frame  of  ref¬ 
erence.  Given  a  flat  satellite  picture,  the  procedure 
for  generating  the  cloud  pictures  is  described.  The 
use  of  shading  to  provide  a  more  effective  three-di¬ 
mensional  picture  is  illustrated.  Then  dimensioning 
of  the  images  is  described,  again  from  the  viewer’s 
frame  of  reference. 

The  visual  and  infrared  image  data  were  col¬ 
lected  by  the  ground  station  at  the  Atmospheric 
Science  Department  at  Colorado  State  University. 
The  computing  and  display  equipment  consisted  of  a 
Digital  Equipment  Corporation  PDP-11/60  mini¬ 
computer  and  a  RAMTEC  9050  Color  Image 
System.  The  RAMTEC  has  a  display  of  512-by-512 
pixels  with  each  pixel  storing  an  8-bit  value.  A  video 
lookup  table  is  used  in  the  RAMTEC  to  define  the 
color  associated  with  each  possible  pixel  value.  A 
13-bit  word  stores  each  color  definition.  The  lower  12 
bits  are  divided  into  3  groups  of  4  bits  each.  Each 
group  controls  the  intensity  of  one  of  the  three  stan¬ 
dard  color  guns  in  the  monitor.  This  provides  a 
possibility  of  16  intensity  levels  for  each  color  gun 
and  a  total  of  4096  colors  and  intensities.  Only  256 
colors  and  intensities  may  be  defined  at  one  time, 
however,  and  only  16  intensities  can  be  generated  for 
a  single  color.  A  more  detailed  description  of  the 
work  is  given  by  Johnson.^ 

Two  Dimensional  Representations 

The  data  that  were  received  from  the  satellite 
consisted  of  visible  and  infrared  data.  In  the  visual 
data,  the  pixel  values  are  related  to  brightness  which 
does  not  directly  relate  to  altitude.  This  limits  the 
usefulness  of  visible  data  for  representing  clouds  in 
three  dimensions.  The  infrared  temperature  data  can 

S3 


Figure  /.  Original  infrared  inwjie. 


be  correlated  to  altitude  via  data  collected  from  bal¬ 
loons  so  that  reasonably  accurate  altitude  represen¬ 
tations  can  be  achieved. 

Thus,  the  first  task  is  to  use  balloon  data  to  re¬ 
late  each  infrared  pixel  value  to  altitude.  After  this  is 
done,  a  new  image  is  generated  to  scale  the  altitude 
values  between  the  values  of  0  and  255  which  repre¬ 
sent  the  lowest  and  highest  altitudes  respectively.  The 
resulting  new  image  can  then  be  displayed  again  as  a 
flat  image. 

An  example  of  the  results  is  shown  in  Figures  I, 
2  and  3.  Given  the  original  image  as  shown  in  Figure 
1;  a  histogram  plot  of  the  data  after  conversion  is 
shown  in  Figure  2.  The  relationship  between  pixel 
values  and  altitude  is  given  by 

A=  (64.5)1  (I) 

where  I  is  the  pixel  value  and  A  is  the  altitude  in 
meters.  The  resulting  reconstructed  flat  image  is 
shown  in  Figure  3  with  a  gray  scale  from  0  to  200. 
Two  hundred  was  used  as  the  upper  bound  since 
values  above  200  are  constant.  A  pseudocolored 
representation  of  Figure  3  is  shown  in  Figure  4.  This 
dramatically  brings  out  the  two  clouds  in  the  center 
of  the  picture  and  these  are  the  two  clouds  we  will 
place  into  a  three-dimensional  image  in  the  next  sec¬ 
tion. 
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To  summarize,  given  a  flat  infrared  image  from 
the  satellite,  a  conversion  to  altitude  is  made  using 
the  best  available  information  to  convert  temperature 
values  to  altitude.  Histogram  values  are  then  ad¬ 
justed  between  0  and  255.  This  may  cause  distortion, 
but  the  impact  of  this  is  still  not  known.  Then  a  new 
image  is  constructed  and  pseudo-colored  if  it  is  to 
show  relevant  features. 
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Hgtire  2.  Pixel  histogram  plot:  pixels  convened  lo  ahitude. 
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Figure  3.  The  image  of  Figure  I  adjusted  for  altitude. 


Three-Dimensional  Representations 

Given  the  infrared  images,  clouds  in  a  desired 
sector  can  be  converted  to  three-dimensional  objects 
that  can  be  viewed  from  any  angle.  The  resulting  im¬ 
ages  can  be  shaded  and  dimensioned  to  provide  the 
viewer  with  a  realistic  view  of  a  cloud  or  group  of 
clouds  from  any  desired  viewing  angle. 

A  surface  can  be  created  from  the  reconstructed 
data  by  means  of  a  collection  of  surface  drawing 
subroutines  that  are  available  from  National  Center 
for  Atmospheric  Research  (NCAR)’  in  Boulder,  Col¬ 
orado  and  were  modified  at  CSU  for  use  with  our 
computer  equipment.  With  the  PDF- 11/60  and 
RAMTEC  system  we  are  limited  to  an  array  of 
32-by-32  points.  Since  the  number  of  pixels  in  the 
portion  of  the  image  used  to  create  the  surface  is 
much  greater,  a  method  to  utilize  the  data  was 
needed.  The  method  which  we  chose  was  to  select 
pixels  at  regular  intervals  so  that  they  form  32  rows 
of  32  columns  within  a  selected  portion  of  the  image. 
The  resulting  surface  can  be  generated  as  a  true  con¬ 
tour  map  or  using  an  X,  Y,  Z  coordinate  system. 
Modification  of  the  programs  resulted  in  surfaces 
with  dimensions.  Here,  labeled  axes  are  used  to  give 
the  dimensions  of  the  X  and  Y  axes  and  a  color  bar 
provides  the  altitude  or  Z  axis  dimensions.  The  soft¬ 


ware  operates  as  follows. 

The  program  is  initiated  by  having  the  user  spec¬ 
ify  the  file  name  of  the  desired  image.  This  image  is 
then  transferred  to  the  RAMTEC,  and  the  lookup 
table  is  defined  for  a  gray  scale  from  8  to  255.  The 
user  interactively  indicates  what  portion  of  this  image 
is  to  be  used  to  create  a  surface.  As  before,  software 
is  used  to  reduce  this  portion  to  an  array  of  32-by-32 
points.  There  are  a  number  of  ways  to  do  this.  Aver¬ 
aging  is  used  in  this  case.  The  portion  of  this  image 
specified  by  the  user  is  selected  so  that  the  number  of 
pixels  in  both  directions  is  a  multiple  of  32.  This  new 
portion  of  the  image  is  then  divided  into  32-by-32 
sections.  The  pixel  values  in  each  section  are  totaled, 
and  these  totals  are  divided  by  the  number  of  pixels 
per  section.  A  potentially  undesirable  effect  of  this 
method  of  reduction  is  low-pass  filtering.  High- 
frequency  features  in  the  image  will  be  attenuated.  It 
should  be  noted,  however,  that  some  information 
will  be  lost  in  any  type  of  reduction. 

The  program  next  converts  the  array  of  infrared 
data  to  altitude  as  discussed  previously.  The  conver¬ 
sion  table  is  created  sometime  before  running  this 
program  and  is  stored  in  a  data  file  which  is  accessed 
to  make  the  conversion  to  altitude.  The  user  specifies 
a  maximum  altitude  that  will  control  the  size  of  the 
cloud  surface  in  the  vertical  direction.  The  distance 
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from  the  zero  elevation  to  the  maximum  altitude  will 
always  be  the  same  for  a  given  orientation. 

Since  ground  features  are  not  necessarily  of  any 
interest  and  tend  to  create  a  more  confusing  surface, 
the  program  now  removes  all  features  below  an  alti¬ 
tude  specified  by  the  user  by  setting  these  array 
elements  to  zero.  The  procedure  up  to  this  point  will 
produce  a  ragged  edge  on  the  bottom  of  the  surface. 
To  make  the  lower  edge  of  the  surface  smooth,  the 
program  changes  the  values  of  the  points  around  the 
edge  of  the  nonzero  surface  to  the  lower  altitude  that 
was  specified  by  the  user. 

The  program  then  creates  the  surface  on  the 
RAMTEC  screen  through  the  use  of  display  sub¬ 
routines.  The  surface  draws  the  X,  Y,  and  Z  direc¬ 
tions.  User  interaction  specifies  how  shape-defining 
sets  of  contours  are  to  be  drawn.  Each  set  of  con¬ 
tours  is  chosen  to  be  close  together  and  cover  a  nar¬ 
row  band  of  altitudes.  The  surface  is  drawn  over  and 
over  again  each  time  with  a  different  set  of  contours. 
Each  set  of  contours  is  drawn  with  a  different  pixel 
value  which  is  defined  as  a  different  color.  The  result 
in  this  program  is  a  surface  on  which  the  color 
changes  from  blue  at  low  altitudes  to  white  at  high 
altitudes.  This  produces  a  known  relationship  be¬ 
tween  color  and  altitude  within  the  program.  The  en¬ 
tire  surface,  including  points  that  are  zero,  is  drawn 
one  last  time  without  any  contours.  The  pixel  value 
used  is  defined  as  black.  The  corners  are  made  equal 


to  zero  to  provide  screen  coordinates  that  will  be  used 
to  draw  the  axes.  The  angle  of  view  for  the  surface  is 
specified  by  the  user. 

Last  of  all,  this  program  adds  dimensions  to  the 
surface  if  so  desired.  The  horizontal  directions  will  be 
described  with  labeled  axes.  Vertical  size  will  be  in¬ 
dicated  by  a  color  bar  with  altitudes  marked  on  it. 
The  corners  of  the  surface  that  were  set  to  zero  pro¬ 
vide  the  screen  coordinates  for  the  ends  of  the  axes. 
The  user  must  then  interactively  choose  three  of  the 
four  corners  to  define  the  end-points  of  the  horizon¬ 
tal  axes.  The  first  corner  chosen  is  the  one  with  a 
screen  coordinate  that  is  lowest  on  the  screen.  The 
other  two  corners  are  those  adjacent  to  the  first. 
These  are  identified  as  right  and  left  by  comparing 
the  X  screen  coordinates.  Two  lines  are  draw  n  on  the 
screen  for  the  axes  using  these  screen  coordinates.  In 
order  to  find  the  scaling  for  each  axis  it  must  be 
determined  whether  the  axis  represents  north-south 
or  east-west.  In  this  case,  if  the  subscripts  of  the  first 
corner  are  equal,  then  the  axis  on  the  right  represents 
an  east-west  direction. 

The  distance  represented  by  an  axis  can  be  found 
by  multiplying  the  distance  per  pixel  by  the  number 
of  pixels  used  to  create  the  surface  in  that  direction. 
The  distance  per  pixel  must  be  entered  by  the  user 
and  is  dependent  on  the  creation  of  the  original  im¬ 
age  data  file.  Once  the  length  represented  by  an  axis 
is  known,  the  interval  between  tick  marks  can  be 
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chosen.  This  program  chooses  ihe  interval  to  be  some 
power  of  ten.  This  power  is  found  by 

p  -  ini(log,„(l.))  (2) 

where  lOf  is  the  interval  represented  between  tick 
marks  and  1.  is  the  total  length  represented  by  the 
axis.  The  distance  between  tick  marks  in  terms  of 
screen  coordinates  is  found  by 


where  s  is  the  distance  between  tick  marks  in  screen 
coordinates  in  the  .X  or  Y  direction,  A  is  the  total 
length  of  the  axis  in  screen  -oordinates  in  the  .X  or  Y 
direction,  T  is  the  distance  represented  by  the  tick 
mark  interval  and  L  is  the  total  distance  represented 
by  the  axis.  By  knowing  the  tick  mark  interval  in 
screen  coordinates,  one  can  use  a  loop  to  draw  tick 
marks  at  these  intervals  until  a  location  is  found  for  a 
tick  mark  that  is  not  on  the  axis.  The  values  written 
on  the  screen  next  to  these  tick  marks  are  multiples  of 
the  interval.  The  value  at  the  first  corner  is  zero.  This 
procedure  is  done  for  each  axis. 

To  indicate  vertical  dimensions,  a  vertical  color 
bar  is  drawn  in  the  upper-right  corner.  The  color  on 
this  bar  changes  from  bottom  to  top  in  the  same  way 
that  the  surface  does.  The  relationship  between  the 
pixel  value  u.sed  to  create  a  color  on  the  color  bar  and 
its  representative  altitude  was  stored  when  the  sur¬ 
face  itself  was  drawn.  The  altitudes  at  the  color 
boundaries  are  written  on  the  screen.  The  program 
then  writes  the  units  of  measurement  on  the  screen. 

Examples  of  the  results  are  shown  in  Figures  5 
and  6  which  are  three  dimensional  views  of  the  two 
clouds  show  n  in  the  center  of  Figure  4.  Figure  5  is  an 
intermediate  result  with  the  flat  base  still  attached. 
Figure  6  is  the  finished  product  with  some  shading 
employed  and  dimensions  included.  Note  that  the  im¬ 
age  in  Figure  6  is  viewed  from  a  different  position 
than  the  image  in  Figure  5,  demonstrating  the  option 
of  different  viewing  angles. 

A  second  method  for  viewing  clouds  as  three- 
dimensional  objects  utilizes  a  set  of  polygons  to  ap¬ 
proximate  the  three-dimensional  surface.  A  software 
package  called  MOVIE,  developed  at  Brigham 
Voung  University,  is  available  as  a  starter  for  doing 
this.'  The  MOVIE  package  contains  the  initial  sub¬ 
routines  needed  for  polygon  approximation  of  a  sur¬ 
face.  arbitrary  coloring  and  shading  of  the  resulting 
surface,  and  the  use  of  transparency  to  allow  the 
viewer  to  see  objects  in  back  of  other  objects.  At 
present.  MOVIE  does  not  run  well  on  the  16  bit 
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/inure  .5.  /niermediale  three-dimensional  represenaiion  of  the 
two  clouds  in  center  of  /inure  4. 


PDP-1 1/60  computer  because  of  the  lack  of  dynamic 
memory.  However,  we  have  been  able  to  obtain  sim¬ 
plified  approximations  to  demonstrate  the  use  of  the 
software  for  cloud  representations.  The  dimension¬ 
ing  software  was  programmed  for  inclusion  in 
MOVIE. 

As  before,  the  initial  image  must  be  sectored  and 
operated  on  to  obtain  altitude  information.  Because 
of  the  lack  of  memory  only  crude  polygon  approx¬ 
imations  can  be  generated,  but  these  are  often  suffi¬ 
cient  for  cloud  representations.  An  example  of  the 
results  is  shown  for  the  previous  two  clouds  in 
Figure  7.  We  call  these  slab  cloud  representations. 

W'hile  the  polygon  approximations  appear  to 
give  a  much  coarser  representation,  the  use  of  a 
larger  computer  will  allow  a  much  better  polygon  ap¬ 
proximation  to  be  obtained.  We  anticipate  'he  results 
will  be  more  cloud-like  although  this  hypothesis  must 
be  proven. 

As  a  final  example,  we  have  recently  used  our 
software  to  depict  the  output  of  computer  models  for 
modeling  cloud  growth.  Three  of  the  results  are 
shown  in  Figures  8,  9  and  10  which  show  the  begin¬ 
ning,  an  intermediate  stage  and  the  final  stage  of  a 
simulation 


Figure  6.  Three-dimensional  cloud  picture  with  dimensions. 
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Conclusions  and  Future  Work 


This  paper  lieseribes  il'.e  methods  fur  generating 
three-dimensional  eloud  images  from  data  generated 
by  geosynehronous  weather  satellites.  The  resulting 
images  are  reasonable  representations  of  elouds  from 
an  observer’s  point  of  view  .  The  use  of  X,  Y,  and  /. 
dimensions  prot  ides  the  viewer  with  accurate  inl'or- 
mation  coneerning  eloud  si/e  from  the  observer  view¬ 
ing  position.  The  eloud  images  ean  be  easily  manipu¬ 
lated  to  permit  observation  from  any  desired  viewing 
angle  and  from  any  distanee.  T  he  use  of  imaging  to 
represent  imxlel  output  is  also  deseribed. 

The  work  was  done  using  a  small  computer  sys¬ 
tem.  The  next  phase  will  utilize  a  DT.C'  VAX-1 1/780 
computer  and  this  w  ill  greatly  enhance  our  ability  to 
investigate  more  sophisticated  ways  of  utilizing  and 
presenting  three-dimensional  data.  By  utilizing  the 
full  version  of  MOVIE,  we  w  ill  be  able  to  more  effec¬ 
tively  represent  clouds  using  better  polygon  approx¬ 
imations  to  solids  and  improved  shading  algorithms. 
The  use  of  a  transparency  algorithm  will  allow  us  to 
"see  through”  clouds  and  to  study  cloud-cloud  inter¬ 
actions.  The  transparency  capability  will  allow  the 
modeling  personnel  to  see  the  inside  of  their  model 
generated  clouds,  for  example,  stream  lines  can  be 
displayed. 

Work  will  also  be  started  to  relate  cloud  data  to 
model  output.  This  will  be  done  by  fTst  displaying 
real  cToud.s  and  model  clouds  in  three  dimension.,. 
From  this  we  hope  to  learn  how  to  use  satellite  data 
better  for  model  initialization.  We  will  also  begin  to 
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/  i^urc  Cloud  represvniuiion  usina  a  crude  polygon  up 
proxtmutton. 

use  temporal  three-dimensional  cloud  images  to 
study  methods  for  estimating  both  horizontal  and 
vertical  wind  fields.  These  studies  will  utilize  our 
present  man-interactive  computing  facilities  to  allow 
our  scientists  to  evaluate  quickly  algorithmic  ap¬ 
proaches  to  the  problems.  ■ 
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Atmospheric  Boundary-Layer 
Processes  Over  Coastal  Waters 
and  Across  Oceanic  Fronts 

Supported  mainly  by  the  Coastal  Sciences  Pro¬ 
gram  of  ONR  and  Naval  Environmental  Predic¬ 
tion  Research  Facility,  Dr.  S.  A.  Hsu  and  his  col¬ 
leagues  at  the  Coastal  Studies  Institute  of  Louisiana 
State  University  are  studying  the  atmospheric 
changes  across  oceanic  fronts  as  well  as  mesoscale 
processes  associated  with  airflow  adjustment  over 
coastal  waters.  When  cold  air  is  situated  over  warm 
water  the  convective  velocity  scale  and  the  height  of 
the  mixing  layer  can  be  estimated  by  the  surface  heat 
flux.  On  the  other  hand,  when  warm  air  flows  over  a 
cooler  sea,  the  height  of  the  modified  layer  is  propor¬ 
tional  to  the  square  root  of  the  fetch  downwind  from 
the  coast.  It  can  be  shown  that  this  equation  is 
superior  to  that  based  on  dimensional  analysis.  All  of 
the  above  concepts  have  been  verified  by  field  experi¬ 
ments.  At  times,  branches  of  warm  oceanic  currents 
such  as  the  Gulf  Stream  and  the  Kuroshio  Current,  as 
well  as  cold  oceanic  eddies  associated  with  upwelling 
such  as  the  Somali  Current,  are  very  near  the  coast. 
Their  impact  on  the  atmospheric  boundary  layer  is 
studied.  It  has  been  demonstrated  that  there  is  ap¬ 
proximately  a  30<t'o  increase  in  relative  humidity  and 
mixing  ratio  from  the  cold  side  to  the  warm  side 
about  15km  downwind  from  the  oceanic  front  in  the 
Korea  Strait.  Also,  a  140m  increase  in  the  mixing 
height  from  the  cold  to  the  warm  side  has  been 
observed.  An  example  is  an  inversion  layer  that  ex¬ 
isted  between  approximately  200m  and  600m  over  a 
cool  pool  because  of  the  presence  of  a  cold  eddy  off¬ 
shore  of  the  Somalia  coast  in  summer;  no  inversion 
layer  was  observed  on  the  warm  side  of  the  upwelling 
zone.  ■ 

(Dennis  Conlon,  ONR) 


Development  of  Techniques  for  the  Remote 
Sensing  of  Surface  Currents  and  Current 
Shear  using  Microwave  Radar 

Surface  currents  are  defined  as  the  mean  hori¬ 
zontal  flow  of  water  within  a  layer  which  extends 
downward  from  the  surface  to  a  depth  of  several 
meters.  The  accuracy  of  conventional  moored  cur¬ 
rent  meter  measurements  suffers  in  this  layer  because 
of  competing  wave  amplitude  and  velocity  effects. 
Lagrangian  drifter  measurements  of  spar  buoys  pre¬ 
sent  difficult  tracking  and  deployment  problems  and 
dye  studies  of  surface  flow  are  difficult  to  interpret 
quantitatively.  Still  the  need  persists  to  accurately 
map  surface  currents,  because  anything  that  floats  on 
the  surface  (ships,  pollutants,  icebergs,  nutrients  and 
even  waves)  is  transported  by  these  currents. 

Recently,  a  dual-frequency  microwave  radar 
technique  has  been  developed  at  the  Naval  Research 
Laboratory  which  is  capable  of  remotely  sensing 
both  ocean  surface  currents  and  vertical  profiles  of 
current  strength  to  a  depth  of  several  meters.  The 
radar  measures  currents  by  monitoring  the  perturba¬ 
tions  in  ocean  wave  propagation  velocity  caused  by 
the  presence  of  the  current.  This  sensitive  measure¬ 
ment  is  made  possible  by  the  dual-frequency  radar’s 
unique  ability  to  “resonate”  with  one  ocean  wave  at 
a  time.  The  particular  wave  that  is  being  used  as  a 
current  measurement-tool  can  be  changed  at  will  be 
changing  the  separation  between  the  pair  of  radar 
frequencies  which  are  scattered  from  the  ocean’s  sur¬ 
face.  Temporal  and  spatial  (mapping)  studies  of  sur¬ 
face  currents  have  been  carried  out  using  this  basic 
mode  of  operation.  Modern  microwave  technology 
facilitates  an  expansion  of  the  radar’s  capability.  In 
particular,  signal  multiplexing  techniques  using 
PROMS  (programmable  read-only-memories)  have 
been  developed  which  enable  the  radar  to  examine  si¬ 
multaneously  the  effect  of  currents  on  a  set  of  up  to 
one-hundred  independent  ocean  waves.  Because  long 
wavelength  waves  are  affected  by  much  deeper  cur¬ 
rents  than  shorter  wavelength  waves,  the  multiplexed 
data  can  be  analyzed  to  obtain  information  on  the 
vertical  distribution  of  currents  beneath  the  ocean’s 
surface.  Information  about  currents  near  the  surface 
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